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We develop a theory for nucleation on top of islands in epitaxial growth based on the derivation 
, of lifetimes and rates governing individual microscopic processes. These in particular include the 

■ encounter rate of j atoms in a state, where in total n>j atoms are present on top of the island, 
' and the hfetime of this state. The latter depends strongly on the additional step edge barrier A_Es 

for descending atoms. We present two analytical approaches complemented by kinetic Monte Carlo 

■ simulations. In the first approach we employ a simplified stochastic description, which allows us to 
'~i I derive the nucleation rate on top of islands explicitly, if the dissociation times of unstable clusters 

. can be neglected. We find that for small critical nuclei of size i <2 the nucleation is governed by 

' fluctuations, during which by chance i+1 atoms are present on the island. For large critical nuclei 

i > 3 by contrast, the nucleation process can be described in a mean-field type manner, which for large 
' ^ ^ ' ■ AEs corresponds to the approach developed by Tersoff et al. [Phys. Rev. Lett. 72, 266 (1994)]. In 

' both the fluctuation-dominated and the mean-field case, various scaling regimes are identified, where 

, the typical island size at the onset of nucleation shows a power law in dependence on the adatom 

1— ( ' diffusion rates, the incoming atom flux, and the step edge crossing probability exp{—AEs/kBT). 

' Although it is possible to extend the simplified approach to more general situations, its applicability 

^ . is limited, if dissociation rates of metastable clusters enter the problem as additional parameters. 

• ' For such situations the second semi-analytical approach becomes superior. This approach is based 

, on novel rate equations , which can easily be solved numerically. Both theoretical approaches yield 

. good agreement with Monte Carlo data. Implications for various applications are pointed out. 



a recent review on both kinetically and thermodynami- 
cally induced instabilities in MBE see ref. During the 
MBE experiment, two-dimensional islands composed of 
adsorbate atoms form on the substrate. If these islands 
coalesce before stable clusters nucleate on top of the is- 
lands in the second layer, a flat two-dimensional film re- 
sults. By contrast, if the onset of second layer nucleation 
precedes island coalescence, then three-dimensional clus- 
ter formation is obtained. The term "second layer nu- 
cleation" should not be taken literally here but rather 
should apply to the formation of stable nuclei composed 
of (i + 1) atoms on top of islands in general. As for the 
equilibrium structures, it might be possible that cluster 
formation sets in above a certain film thickness, when 
the relevant parameters governing the nucleation of sta- 
ble clusters on top of islands (see below) depend sensi- 
tively on the film thickness. However, despite this simi- 
larity of the possible growth processes with the equilib- 
rium growth modes, it should be noted that the dynamic 
problem is very different. In MBE flat films can be pro^ 
duced even if the adsorbate does not wet the substrate.B 

The first theory for second layer nucleation in MBE 
was set up by 1. Tersoff, A. W. Denier van der Gon, and 
R. M. Tromp,LJ which will be referred to as "TDT ap- 
proach" in the following. Solving the stationary diffusion 
equation in the presence of an iacoming flux and employ- 
ing classical nucleation theory,EJ these authors succeeded 



. I. INTRODUCTION 

o 

A fmidamental problem in the theory of thin film 
growtffl is the question, under which condition flat, two- 
dimensional films form on the substrate surface in con- 
^ trast to mutually-separated, three-dimensional clusters. 

, For films growing under equilibrium conditions, this ques- 
^ID tion has been answered many years ago: If the inter- 
_ facial tension between the substrate and adsorbate is 
. larger than the difference of the respective surface free 

■ energies, then duster formation is preferred ("Volmer- 

■ Weber growth" 0), while a smaller (or equal) interfacial 
tension leads to the formation of flat films ("Van-der 
Merwe growth"). An inteonediate case is the "Stranski- 

^ ■ Krastanov growth" mode,0 where cluster formation sets 
' ■ in after the thickness of an initially smooth film exceeds a 
^ critical height. This case may be understood from an in- 
Q terfacial tension that varies with the film thickness. More 
^ : recently, the influence of strain effects on equilibrium fil: 
morphologies has been investigated by various authors 
Films developing in the process of molecular beam 
' epitaxy (MBE) are usually not in thermodynamic 
equilibrium.^. Rapid growth of films is achieved by a high 
supersaturation of the vapor, and growth kinetics is gov- 
erned by evaporation, diffusion and aggregation processes 
far from equilibrium. Determining the film morphology 
in this situation is a problem of stochastic dynamics. (For 



* Current address: Department of Physics and Astronomy, The Johns Hopkins University, 
3400 N. Charles Street, Baltimore, MD 21218, U.S.A. 



1 



in deriving an explicit expression for the rate of nucle- 
ation ^l{R) on top of circular shaped islands of radius R. 
Assuming all island radii to evolve approximately as the 
mean island radius R{t) ~ t^^^ at time t (this situation 
will be referred to as the "single-island model" in the 
following), they calculated the fraction fo{t) of "covered 
islands" (i.e. on top of which a stable cluster has nucle- 
ated) from ri(i?). It turned out that fo{t) rises from zero 
to one in the vicinity of a "critical time" tc, which allows 
one to define a critical island radius Rc = R{tc) for second 
layer nucleation. A simple criterion for the occurrence of 
"rough multilayer" as opposed to smooth "layer-by-layer 
growth" is that Rc is smaller than the mean distance I 
between islands in the first layer. 

An important factor controlling the film morphology is 
the additional step, edge barrier AEs—Es — Eq (Ehrlich- 
Schwoebel barriem) that has to be surmounted by an 
adatom in addition to the bare surface diffusion barrier 
Eq, when an adatom crosses an island edge.EI For larger 
AEs one expects adatoms to remain longer on islands 
and therefore to accumulate more easily, which would 
lead to an increased second layer nucleation rate ^1{R) 
and a smaller Rc- In fact, the theory predicts that only 
for sufhciently large AEs three-dimensional clusters can 
occur on the substrate (for an alternative possibility see 
however ref . |l^) . Using the TDT approach,, rAEs was 
estimated for a variety of different systems .llirL3 

An alternative approach for treating the problem of 
second layer nucleation within a stochastic description 
based on scaling arguments was developed recently by 
us.li3 It was shown that for i = 1 the TDT approach 
is not applicable, but the detailed treatment of fluctua- 
tions with only two atoms on top of the island yields a 
correct description of the process (see also ref. ^Tf for an 
earlier approach focusing on one dimension see ref. [l8[ ). 
In this work we will extend our former study of sec- 
ond layer nucleation by means of both kinetic Monte 
Carlo simulations and scaling analysis. In particular, 
we will show that the mean-field assumptions underly- 
ing the TDT approach are valid for large critical nuclei 
i > 3, while for small critical nuclei i = 1, 2 second layer 
nucleation is dominated by fluctuations. Furthermore, 
we develop a novel rate equation approach, which allows 
one to calculate the time- development of cluster configu- 
rations on compact two-dimensional islands under quite 
general conditions. 

From the outset, one should distinguish between nucle- 
ation on top of islands with compact shape as opposed to 
nucleation on islands with strongly ramified shape. In the 
latter situation, diffusion of adatoms on the islands be- 
comes a rather complex phenomenon due ta_the confined 
motion along branches of various lengths .Ej We restrict 
our discussion to nucleation on compact islands here. In- 
tuitively, one would expect second layer nucleation on 
ramified islands to be unlikely, so that this restriction is 
of minor importance. Moreover, it should be noted that 
even for compact shapes, the island bowndaries may have 
a fractal or, more precisely, self-afRnec2l structure (this is 



the case e.g. for Eden clusterB). Then the microscopic 
step edge barrieiE3 can vary strongly along the island 
boundary. In any case we will always understand AEs 
as an effective barrier (see below and ref. |2^) . 

Parameters governing the second layer nucleation are 
the incoming atom flux _F, the jump rate D/a? of 
adatoms, the step edge barrier AEs, and various dissoci- 
ation rates of unstable clusters of size s<i. If the bond 
energies of the unstable clusters (i.e. of clusters of size 
s<i) are negligibly small, then the nucleation rate Q,{R) 
and critical radius Rc depend only on two dimensionless 
parameters. These are the ratio 
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(with a being the lattice spacing in the substrate plane) 
and the edge crossing probability 



a = exp 



A^s 
fceT 
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Based on a simplified stochastic description (see also 
ref. |l|) we will argue that for small critical nuclei i = 1, 2 
the mean number of atoms on top of the island is smaller 
than one and the stable nucleus is formed due to fluctu- 
ations. This gives rise to four scaling regimes in an a — F 
diagram, where Rc ^ V^a^ with different exponents 7 
and fi. For i > 3 by contrast, nucleation starts out from a 
situation with many atoms present on the island. Under 
these circumstances, three different scaling regimes can 
be identified, and two of them correspond to the ones 
predicted by the TDT approach. By comparing Rc with 
the mean distance I between islands on the substrate sur- 
face, the transition line separating rough multilayer from 
smooth layer-by-layer growth is identifled in the a — T 
diagram. When the bond energies of unstable clusters 
become appreciable, the corresponding dissociation rates 
enter the problem as additional relevant parameters. It 
becomes difficult then to separate scaling regimes in prac- 
tice, and the simplified stochastic description becomes 
of limited value. However, by employing the novel rate 
equation approach it is still possible to determine /o(t) 
and Rc in a simple manner. 

Moreover, we will discuss how to derive the fraction of 
f(t) of covered islands, when one relaxes the assumption 
that all island radii evolve as the mean radius R{t). In 
the time regime of almost constant island density ( "satu- 
ration regime" preceeding island coalescencecj) we can 
define an effective "capture_area" for adatoms by the 
Voronoi cell for each island.LJ In order to calculate f{t) 
for the "multi-island model" from fo{t), one needs to 
know the probability distribution of islands with a cer- 
tain size and capture area, when the saturation regime is 
reached. 

The paper is organized as follows. In Sec. |ll| we give 
a short review on basic concepts used in the descrip- 
tion of submonolayer growth and discuss in part II A im- 



portant quantities and equations underlying the physical 
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processes involved. In part II B we summarise the results 



of the TDT approach. We then continue with a detailed 
description of the simulation techniques in Sec. [II and 



test the equivalence of the multi-island and single-island 
model. 



In Sec. IV a simplified stochastic description of second 



layer nucleation is first presented in its general methodol- 
ogy and subsequently employed to small and large critical 
nuclei. It is instructive to observe the different physical 
conditions that lead to the nucleation event in the two 
cases. After showing that the treatment of metastable 
clusters is rather complicated within the simplified de- 
scription, we discuss in Sec. ^ a general approach for 
second layer nucleation on the basis of novel rate equa- 
tions. Sec. VI concludes the paper with a summary and 



discussion of the most important results as well as an 
outlook to further research. 



II. BASIC QUANTITIES AND CONCEPTS 
A. Atomistic processes in thin film growth 

In MBE atoms are deposited on a substrate surface 
with a rate Fa^ per unit cell. At very high tempera- 
tures the adatoms reevaporate but under ordinary con- 
ditions this reevaporation can be neglected or effectively 
taken into account by a reduced deposition rate. Once 
an adatom is deposited, it starts a thermally activated 
diffusive motion with jump rate D/o? oc exp(— .Eq/'sbT')- 
Adatoms come into contact as time progresses, and form 
islands that are held together by some bonding energy. 
While unstable islands of small size s < i dissociate, is- 
lands with size s > i are stable (on all relevant time scales 
of the experiment). An island of size s = i is called a crit- 
ical nucleus .0 

Islands of larger size are formed by aggregation of 
adatoms (or small mobile islands) to existing immobile 
islands, and by coalescence. These processes lead to an 
island size distribution which becomes broader with in- 
creasing time. In the submonolayer regime, the most 
important physical questions are: (i) How large is the 
density px {t) of stable islands on the substrate surface at 
time tl (a) What is the form of the distribution Xs(s,i) 
of island sizes s at time t (with s being the number of 
atoms forming the island)? (iii) What do the stable is- 
lands look like? These questions have been extensively 
studied in the past, both by experiment and by theory. 
We will briefly summarise those results, which are rele- 
vant for the following analysis. 

The typical behavior of Px{i) is depicted in Fig. |^. 
Also shown is the ad|atom density pi [t] . As suggested by 
Amar and FamilyJHJ one may distinguish between four 
different time regimes: The low-coverage regime L where 
Px{t) increases with t and Px{t) < Pi{t), the intermedi- 
ate coverage regime I, where Px{t) increases with t and 
Px{t) > Pi{t), the saturation regime S (called aggregation 
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FIG. 1. The densities pi{t) and px{t) of adatoms and sta- 
ble islands as obtained in a kinetic Monte Carlo simulation 
modeling submonolayer growth (for an analogous figure see 
e.g. ref. M). The different time regimes of low coverage L, 
intermediate coverage I, saturation S and coalescence C are 
indicated. 

regime A in ref. , where Px {t) stays approximately con- 
stant, and the coalescence regime C, where Pxit) strongly 
decreases due to coalescence of stable clusters. At high 
coverages 9 = Fa^t, the monomer density pi{t) becomes 
small (see Fig. 1), and the standard rate equations for 
submonolayer growthEl predict px to evolve as 

_ i 

Px{t)oc(-^y'^\Fah)The^^/(^+'^'^-^ 



(3) 



where Ei is the bonding energy of the critical nucleus in 
its preferred atomic configuration. 

It should be noted that eq. (||) is not valid in the satu- 
ration regime S, where the densities of islands with sub- 
critical and critical size are very small (unless there are 
metastable subcritical nuclei). In this regime almost all 
adatoms being deposited attach to preexisting stable is- 
lands, so that Px{t) stays constant, Px{t) = Px- Within 
the standard rate equation approach, this effect may be 
accounted for by a proper dependence of the "capture 
numbers" on the adatom density pi{t) (for a detailed 
discussion of this point in relation to experiments see 
ref. The scaling of px with F, however, is still cor- 
rect in regime S, 



(4) 



One of the most detailed studies of the island siae 
distribution has been perfopned by Amar and FamilyEfl 
based on the scaling ansatznS^Ej 



Xs{s,t) 



■f 
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Here Xsis,t) ~ ps{t) / ptot{t) is the probability that an 
island has size s; (psit) is the density of islands with size 
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edge diffusion is suppressed, one obtains dendritic or ran- 
dom fractal structures. Dendritic growth is preferred 
at low T or small F, and a shape transition from den- 



dritic to random fractal structures has been fouK 



e.g. for 

An ex- 



Ag/Pt(lll) upon lowering the deposition flux 
ample for a random fractal structure obtained in a com- 
puter simulation is shown in Fig. ^ja. At high temper- 
atures, edge diffusion becomes relevant, and polygonal 
or "irregular" compact island morphologies develop (see 
e.g. ref. |3^ ). An example for an irregular structure is 
shown in Fig. |^. 

As mentioned in the Introduction, for the study of 
second layer nucleation we will focus on compact island 
shapes. Moreover, second layer nucleation in the inter- 
mediate regime I is unlikely to occur, since the island 
radii in this regime are typically smaller than the crit- 
ical radius i?c- We therefore consider the second layer 
nucleation in the saturation regime S, where eqs. (^j^) 
apply. 



FIG. 2. Typical island morphologies as obtained in kinetic 
Monte Carlo simulations (a) for strictly irreversible attach- 
ment and (b) for irreversible attachment with local relaxation, 
in the absence of preferred growth directions. For the simula- 
tion technique of the local relaxation process see Sec. ^ and 
Fig. I 

s and ptot{t) is the total island density). (. . .) denotes an 
average over s with respect to Xsis,t). From (0) fol- 
lows, for (s™) < oo, (s™) = (s)"" dxx'^fix), and 
taking m = 0, 1 one obtains J dxf{x) — Jdxxf{x) = 
1. The mean island size {s{t)) is given by (s(i)) = 
Y^7LiSPs{t)/ptotit) = 6 / ptot{t)a'^ ■ In the saturation 
regime S, in particular, where ptot — Px, since the den- 
sity of islands with subcritical and critical size is small, 
it follows from eq. 



{sit)) 



Pxa" 



)F^ 



(6) 



The relation {sit)) ~ theta/ pxO? — Ftj p^ can also be un- 
derstood more directly, since the increase of (s(i)) with 
t is given by the flux times the mean capture area of 
adatoms. 

The, scaling function /(m) was suggested to have the 
forirO 



/(u) = CiU*exp(-WiU^/''") 



(7) 



in regime S. This function has a maximum at u = 1 and 
the two conditions J dxf{x) = J dxxf{x) ~ 1 determine 
the parameters Ci and a^. Equations (^[|^) have been 
shown to give a fairly good approximation of some sim- 
ulations and experiments. 

The protJjCm of the island shapes is not yet well 
understood,Q but one may roughly answer the third ques- 
tion posed above as follows. For strictly irreversible at- 
tachment, where local relaxation of atoms due to fast 



B. TDT Approach 

In the TDT approach,! one starts by calculating the 
adatom density pf on a circular island with radius R 
in the stationary state. The stationary diffusion equa- 
tion with the incoming atom flux acting as a source term 
reads 
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and it is supplemented by the boundary conditions [a 
expi-AEs/kBT)) 



dpf 
dr 



r=0 



0, 



dpf 
dr 



a 



r=R 



r=R 



(9) 



where aJ-oi is commonly referred to as the "Schwoebel 
length" .eJ The boundary conditions express the fact that 
the current density —Ddpf/dr must vanish at the origin 
and that at the edge it is given by the density pi times 
the "velocity" (rate times lattice spacing) {Da/a?)a to 
cross the step edge barrier. The solution of eqs. (0,0) is 



Pt{r) = pm^ 



Fr^ 
4D 



, FR^ f 2a 
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According to standard rate equation theory! the local nu- 
cleation rate is proportional to Dp^^ , so we obtain from 
eq. (0) 



Vt{R) = 
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a > 



2a 

'r 



where k is a constant. 

For a given time evolution of the island radius R = R{t), 
one can calculate the probability /o(i) for a stable nu- 
cleus to have formed on top of the island up to time t as 
follows: The increase fo{t + At) — fo{t) in a small time 
interval At is equal to the probability [1 — fo{t)] that up 
to time t no stable nucleus has formed times the prob- 
ability U{R{t))At that the nuclcation takes place in the 
time interval [t,t + At]. Taking the limit At ^ and 
solving the corresponding differential equation with the 
initial condition /o(0) = yields 



fo{t) = 1 - exp 



dt'n{R{t')) 



(12) 



For compact island growth during an MBE experi- 
ment, we have R{t)^ {s(t))^^'^ and thus from eq. (||) 



A(Fa2t)i/2r/2(-*-2) 



(13) 



with A being some constant. Inserting this growth law 
into dri|,pl) yields 



fo{t) = 1 - exp 



2r"^ 



R{t) 



dr rfl(r) 



(14) 



Equations ( |15[jl7| ) predict that for large step edge bar- 
riers, Rc depends strongly on AEs, Rc ~ exp[— (i -|- 
l)AEs/{i + 5)kBT], while for small barriers, Rc becomes 
independent of AEg. For j = 1 in particular, one finds 
Rc ~ r2/9ai/3 for a < T^/^ and i?. ~ T^^^ for a » T^/^. 



III. KINETIC MONTE CARLO SIMULATIONS 



Kinetic Monte Carlo simulations are a weU^steblished 
technique for modeling MBE experiments. In our 
investigation of second layer nucleation, we adopt a sim- 
ulation scheme similar tojpsevious, successful models of 
surface growth kineticsHJ'cS We chose a substrate with 
fcc(lll) symmetry, since surfaces of that kind are of- 
ten studied in metal epitaxy, and commonly exhibit high 
AEs- In a full simulation scheme of the growth kinetics 
(multi- island model), we include all processes of evapo- 
ration, diffusion and aggregation occurring in the MBE 
experiment. By analyzing the set of islands of various 
size on the substrate, we determine the fraction f{t) of 
covered islands at time t. On the other hand, we consider, 
as in the TDT approach, only one island with the mean 
radius R(t) evolving deterministically in time. The frac- 
tion fo{t) of covered islands in this single-island model is 
then determined by calculating the probability for second 
layer nucleation up to time t from a large set of indepen- 
dent simulations. By examining both models we are able 
to quantify the influence of the cluster size distribution 
under generic growth conditions. 



1 — exp 



1 — exp 
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where C> ee {2-^('+^^TrKA-^)/[{i + 2){i + 3)] and C< = 
2^VkA^^. In going from the first to the second line in 
( |l^ ) we have used that the integral over r is dominated 
by the upper integration bound R{t) (for R{t)/a ^ 1). 
It follows that the critical radius scales as 



where 



Rc ^ r^a^ , 



{i + 2){i + b) 
i 



and 



2{i + 2) ' 

(';+i) 



a « r/[2(^2)] 
a » r'/[2(^2)] 



(15) 
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M = 



(« + 5) 

, a » r*/P(H-2)] 
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A. Multi-island model 



Atoms are randomly deposited with a rate Fa^ per 
unit cell onto a triangular lattice. After instantaneously 
relaxing to a position, where they are supported by three 
nearest neighbors in the layer below ("downward tun- 
neling"), the atoms change their position by performing 
thermally activated jumps to a vacant nearest neighbor 
site in the same layer with a rate D/6a^. Only one atom 
is allowed to occupy a given lattice site. We first consider 
a "non-interacting particle model" where all binding en- 
ergies of subcritical clusters of size s < i are neglected. 
This means that (i + l) adatoms have to encounter each 
other on nearest neighbor sites in order to form a sta- 
ble nucleus. Within the single-island model we later will 
also consider finite binding energies of metastable clus- 
ters, which causes various dissociation rates to enter the 
problem as additional parameters. 

Once a stable cluster of size s>i has formed, adatoms 
can attach to it. Compact island morphologies are known 
to emerge if a fast diffusion process is present along island 
edges. Here we model this process similar as in earlier 
approaches (see e.g. refs. |4 39) by including a local re- 
laxation mechanism. In this method an atom being in 
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FIG. 3. Illustration of the local relaxation process when an 
adatom arrives at an island edge. 

contact with at least one nearest neighbor after a jump, 
is immediately transferred to a nearest neighbor site, if 
it can increase its coordination number. This procedure 
is repeated until the atom can no longer increase its local 
coordination (see Fig. ||). 

Interlayer diffusion of atoms deposited onto islands is 
hindered by the Ehrlich-Schwoebel barrier AEs, which 
reduces the jump rate D/6a^ hv the edge crossing prob- 
ability a — exp(— A£'s/fcBr).c3 For computational con- 
venience, we model the crossing by a two-step process in 
the simulation: First, when an atom passes the boundary, 
it remains in the same layer but moves to a place, where 
it is supported by only two atoms underneath. Then 
the atom immediately drops down to the layer below 
and moves to the nearest "stable" site according to the 
local relaxation mechanism introduced above (for simi- 
lar simulations including AEs, see, e.g. |l|,|l|). We do 
not distinguish between crossing of A and B steps and 
have not attempted to model any more realistic scenar- 
ios, as e.g. collectivaj'ea.rrangements of atoms including 
exchange processes. lJcj This is well justified as long as 
one is iaterested in the influence of an effective Schwoebel 
barrier H3 

In the following, we focus on the case i — 1 first. Typi- 
cal film morphologies resulting from the simulations have 
been shown in Fig. ^b. Note that the boundaries of the 
islands are still rough despite the local relaxation mech- 
anism. The fraction f{t) of covered islands as a function 
of the total coverage Fa^t is shown in Fig. ^ for some rep- 
resentative parameters (full symbols). As expected, f{t) 
first is close to zero, then increases strongly in some time 
interval around a "critical time" tc, and finally saturates 
at one. In the inset of Fig. H we show the dependence 
of the mean island radius R(T) = ((s(t))/7r)^/^a on Fa^t 
during the evaporation. In agreement with eq. (p^), we 
find R{t) = A(i^a2t)i/2ri/6 with A = 0.78. 

To be specific, let us define the critical time tc via the 
condition f{tc) — 1/2, and the corresponding critical is- 
land radius Rc by Rc = R{tc), 

f{Q = 1/2 , Rc = Rite) = ^(Fa^i Ji/2r»/2(^2)^ (^8) 

Plots of Rc as a function of a for various fixed F are shown 
in Fig. H (full symbols). With increasing step edge bar- 
rier, i.e. decreasing a, adatoms on average remain longer 
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FIG. 4. The fractions /(t) (full symbols) and /o(i/l-21) 
(open symbols) of covered island as a function of the to- 
tal coverage FcP't for i = 1, a = 10~^, and three different 
r = 10^ (■, □), 10'' (•, o), and 10^ (A, A). The dashed 
lines fitting /o(f/1.21) were calculated from the theoreti- 
cal predictions for the second layer nucleation rate Q{R{t)) 
(eq. (^il)). The solid lines fitting /(t/1.21) were calculated 
according to eq. ( pl| ) with Xct(o") ~ o"^ (see text). The in- 
set displays the time-dependence of the mean island radius 
R{t) = {{s{t))/n)'^^a (see eqs. (||l|)). 
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a=exp[-AEj/kBT] 

FIG. 5. Dependence of the critical island radius Rc on a 
for i = 1 and various T = 10^ (■, □), lO** (•, o), 10^ (A, A), 
and 10* (T, V). Full symbols refer to the results from the 
multi-island model, while open symbols refer to the results 
obtained from the single-island model [Rc = LIR'^). The 
dashed line marks the onset of layer-by-layer growth, and the 
solid line with negative slope marks the crossover between 
regimes I and II; the solid lines fitting the data in regime II 
have slope 1/7. 

on an island and nucleation of stable dimers occurs at 
smaller island radii. Accordingly, Rc decreases with de- 
creasing a (see "regime IT' in the figure) . For very small 
a, however, the step edge barrier is practically never sur- 
mounted and thus is in effect infinitely high. Therefore, 
Rc becomes independent of a ( "regime I" in Fig. I). The 
crossover between the two regimes is marked by the thick 
solid line. The full symbols in Fig. ^ terminate at the 
dashed line a^(r), which marks the onset of island coa- 
lescence. For a > a^(r), islands in the first layer merge 
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before second layer nucleation takes place and Rc can 
no longer determined from the multi-island model. It 
is important to note that the dependence of Rc on a is 
much weaker than predicted by the TDT approach: The 
solid lines in regime II have slope 1/7 corresponding to a 
power law Rc ~ a^/'' rather than Rc ^ a^/"^ as predicted 
by eqs. ( [Tsl , [TtI ). Moreover, regime I does not occur in 
the TDT approach. 



B. Single-island model 

Second layer nucleation can also be addressed in a 
simpler model, which does not attempt to describe the 
entire growth dynamics, but focuses on the decisive 
factors that determine nucleation in the presence of 
the step-edge barrier. In this model the complicated 
nucleation and diffusion-mediated growth of the two- 
dimensional islands, on which the second layer nucle- 
ation takes place, is replaced by letting the radius of 
one circular island expand deterministically in time as 
i?(t)/a = ^(i^a2t)i/2p/2(*t-2), where A is taken from the 
full simulation of the multi-island model. 

The island is embedded in a substrate area large 
enough to accommodate the island at all relevant times. 
Deposition and diffusion of adatoms take place in the 
same manner as in the multi-island model. Atoms in- 
side the island boundary can escape by overcoming the 
step edge barrier. Those atoms that have surmounted 
the barrier or that have been deposited outside the is- 
land boundary are removed from the lattice. Thus the 
single-island model considers the deposition of random 
walkers within a time-dependent, circular boundary that 
is partially reflecting. Due to its greater simplicity, it 
allows for more specific analysis with a larger parame- 
ter space (there is no restriction due to coalescence of 
distinct islands). 

In the non-interacting particle model, the "critical 
event" is to find (i+l) atoms on neighboring lattice sites. 
Analogous to the multi-island model we can define the 
fraction fo{t) of covered islands up to time t. The frac- 
tion now refers to a set of islands obtained in indepen- 
dent simulation runs. All islands in these runs grow with 
the same deterministic growth law. Results for fo{t) are 
shown in Fig. ^ (open symbols) for the same parameters 
as in the multi-island model. Good agreement with /(t) 
is achieved for small times (corresponding to /(i)<l/2), 
when the time in the single-island model is rescaled by 
a constant factor, i.e. /(t)~/o(t') with t' — t/1.21. The 
factor is a consequence of the idealized circular island 
perimeter in the single-island model. In the multi-island 
model by contrast, the islands are far from being per- 
fectly circular (see Fig. 2|). They have rougher edges with 
more boundary sites, which causes adatoms to escape the 
islands more easily and second layer nucleation to occur 
at later times t = 1.21t'. 

At larger times (corresponding to f{t) > 1/2), however. 



f{t) deviates from fo{t') and these deviations become 
more pronounced for larger F. The reason for this dis- 
crepancy is the presence of islands with size much smaller 
than {s{t)) in the multi-island model. Nucleation of sta- 
ble clusters on top of these islands occurs at a later time, 
which causes f{t) to be smaller than foj t') — / o(t/1.21) 
at large t. In fact, we will show in Sec. Ill C] that this 
effect can be accounted for by considering the probabil- 
ity distribution of islands with a certain size and capture 
area. When Rc approaches the mean distance I, coales- 
cences of larger islands also lead to modifications of f{t) 
for t > tc- 

The critical radius i?^ in the single-island model can be 
defined as in the many island model by i?^ = i?(i^), where 
foit'c) = 1/2. Due to the fact that tc = 1-21^^ we expect 
Rc = l.2l^/'^R'^ = I.IR' Resuhs for l.li?;, as a function 
of a are shown in Fig. (open symbols) for the same pa- 
rameters as in the multi-island model (full symbols). As 
can be seen from the figure, there is almost perfect agree- 
ment between both data sets. Moreover, the data for R'^ 
can be obtained also beyond the dashed line marking the 
onset of layer- by-layer growth. Let us also note that, as 
long as one is interested only in R'^ (or Rc = l.li?^), 
one may obtain it even more simply in the single-island 
model (without calculating fo{t)) by determining the av- 
erage radius of the island at the time of the nucleation 
event. 



' dfojt) 
dt^Rit). 



(19) 



Note that dfo{t)/dt is the probability density of the sec- 
ond layer nucleation times and that the average of R{t) 
with respect to dfQ(t)/dt is approximately equal to R{t'c), 
since dfQ(t)/dt is sharply peaked around t'^. 



C. Equivalence of the single-island and the 
multi-island model 

In order to determine f{t) from /o(t) we define by 
tpis, cr, t) ds da the probability for an island to have a size 
in the interval [s, s-l-iis] and a capture area in the interval 
[ct, (T -|- dcr] at time t, where the capture arjea is given by 
the Voronoi cell associated with an islandE3 

Let us consider fo{t) to be a functional of the growth 
law R{t) only, as it is the case, for example, when one 
approximates the second layer nucleation by a Poisson 
process with a time dependent nucleation rate fl{R{t)). 
Then foit) = Go[R{t)] = 1 - exp[~ dt' n{R{t'))] (see 
eq. [l^). In the saturation regime the growth law for an 



island can be written as ttR (t) 



- Fa{t — t-^), where 



tx is the time when the saturation is reached (see Fig. |l|) 
and Sx is the island size at that time. (We restrict our- 
selves to film-morphologies far from coalescence here, so 
that a can be regarded as time-independent.) With the 
specified growth law, the functional Go[R{t)] can be ex- 
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pressed by a function go = goit; SxjCrjtx), and f{t) is 
calculated via 



OO pQO 



f{t)^ / dsx / daip{sx,cr,ty,) gQ{t;Syi,a,tx) ■ (20) 

Jo JQ 

A detailed investigation of the probability distribution 
ip{s,(T, t) is certainly of interest but beyond the scope of 
the present work. A simple idea would be to neglect 
correlations between the stochastic variables s and tr, 
ip{s, a, t) ~ Xs{s, t)xcr{o', t), and to use previously derived 
scaling forms for the island size distribution Xs(sji) (see 
e.g. refs. |35|, p4 | ) a nd the capture area distribution x<y{'^^ t) 
(see e.g. refs. p9 45). 

Here we will follow a simpler approach. Since for typ- 
ical situations we find both Sx and t-^ to be significantly 
smaller than Sc — t^Rc and tc, respectively, we use the 
growth law s{t) = Fat for an island with capture area 
a in the full simulation. The on-top nucleation prob- 
abilities go{t;a) for islands exhibiting different capture 
areas can then be related by a rescaling of time, i.e. 
go{t;ai) = gQ(crit/a2]<J2)- Moreover, since for film mor- 
phologies far from coalescence (a ^ ai,{T)), x(o', t) is ap- 
proximately independent of time, we have /o(i) =5o(^; 
where a = J daxa-i^jt)^ — Px^ ■ Hence, 

/•oc poo - 

f{t)^ daXa{<j)go{t;<^)= daxa{<y)h[^ty (21) 

In this simplified eq. (|T]) knowledge of the nucleation 
rate Vl{R) is not necessary and f{t) can be directly ob- 
tained from /o(t) when Xct(c) is known. 

Writing X(t(o') — a'^'^h{a /d), where J dxh{x) — 
J dxh(x)x = 1, the transformation ( |2l| ) becomes f[t) = 
dxh{x) fo{xt) . For a random distribution of point is- 
lands, we would have h{x) — exp(— a;). However, since 
there is a depletion zone of adatoms near an island, the 
probability for other islands to nucleate in an area close 
to an existing one is reduced and not exponential. For an 
isolated island, dimensional analysis predicts the exten- 
sion ^ of the depletion zone to be of order {D/FY^^. By 

comparing ^ with the mean distance I ~ px ^^"^ ~ r*/^'^'+^^ 
between islands, we expect that h(x) exhibits a large x 
regime with h{x) ~ exp(— x) only for i > 2. For i — \ 
we thus are satisfied with a simple power law ansatz 
h(x) = Cx^ for X < Xi,, where C and follow from 
the two conditions imposed on h{x), and is a fitting 
parameter. 

To test this ansatz we take /o(i) for r = lO*" from Fig. ^ 
(open symbols or dotted lines) and compare/ (t) as cal- 
culated from eq. (H]) (sohd lines in Fig. k|) with the 
corresponding fit) as obtained in the simulation (full 
symbols in Fig. ^). As can be seen from Fig. ^, for 
F = 10^ and F = 10^ a fairly good account of the differ- 
ences between /o(t) and f{t) can be obtained by choosing 
(f) = 2. However, for F=10^ the theoretical curve under- 
estimates the fraction of covered islands at large times 
(where f{t) > 1/2). Better agreement between theory 



and simulation can only be obtained if one would allow 
(j) to depend on F. Alternatively, we have tried an ansatz 
for h{x) similar to that used by Amar and Family for the 
scaling function characterizing the island size distribu- 
tion (see eq. (^). This ansatz yields comparable results, 
but is also not successful in accounting for the changes 
with F. We finally have to note, that at the time when 
submitting the paper, a more-detailed theoretical account 
for Ms, a, t) was published.cS The use of this finding in 
eq.(En) and the comparison of the resulting f{t) with 
Monte Carlo data will be presented elsewhere. 

Having shown that the single-island and multi-island 
models are essentially equivalent, except for differences 
between f{t) from fo{t) for large times that can be at- 
tributed to the island size distribution, we will focus on 
the single-island model in the remaining part of the pa- 
per. 



IV. SECOND LAYER NUCLEATION IN SIMPLE 
SITUATIONS 



In this Section we develop a stochastic description of 
the nucleation process based on the scaling approach 
for second layer nucleation presented in ref. ^ (see also 
ref. [T^ ). The procedure focuses on the non- interacting 
particle model, although formally it is possible to extend 
scaling concepts also to situations, where the lifetimes 
of unstable clusters become important. This was shown 
by Krug et alS-H a nd is discussed in a more general con- 
text in Sec. |IVE| . The treatm ent o f the non- interacting 
particle model outlined in Sec. IV A already captures the 



salient features of the problem in terms of lifetimes, oc- 
cupation probabilities and encounter rates. We will show 
that there exist two possible mechanisms for the forma- 
tion of a stable cluster: In the first case, there is typically 
no atom on top of the island and a stable cluster is formed 
due to fluctuations, in which by chance i + 1 atoms are 
present on the island. In the second case by contrast, 
there are on average more than i+1 atoms on top of the 
island during the formation of a stable cluster so that the 
nucleation process can be described in a mean-field type 
manner. 

It turns out that the fluctuation-dominated case takes 
place for i < 2, while the mean- field situation occurs for 
i > 3. The TDT approach corresponds to the mean-field 
case with the notable supplement that for very large step 
edge barriers one should deal with the time-dependent 
adatom density pi{r,t) (solution of eqs. (@j9|)) to calcu- 
late the nucleation rate ^{R) from eq.([ll[). In the lan- 
guage of critical phenomena, one may regard i = 2 as 
the upper critical size of the critical nucleus above which 
mean-field theory becomes applicable. We have to note 
that the existence of this upper critical size was not per- 
ceived by us in ref. and accordingly, the extension 
of the scaling arguments for the fluctuation-dominated 
situation to i = 3 was not allowed. 
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In the stochastic formulation presented below we will 
develop many of the necessary ingredients for the general 
treatment of second layer nucleation in the next Sec. 0. 
Moreover, it is discussed under which conditions mean- 
field type expressions ex Dp^^"^ for local nucleation rates 
can be used. 



A. General procedure 

In order to determine a second layer nucleation rate 
ri(i?) we start by considering a time interval At{R), dur- 
ing which R{t) does not change significantly. For exam- 
ple, for the generic growth law ( p^ we may require At{R) 
to correspond to a 10% change of R, which would give 
At{R) = 0.21{R/a)'^ /[A^Fa^r^('+^^], i.e. 



AtiR) - F-ir-V(H-2)^2 



(22) 



The nucleation rate Q{R) is the mean number n^^^{R) of 
nucleation events in time At{R) divided by At{R), 



n{R) 



n„uc(i?) 
At{R) 



(23) 



A nucleation event occurs, if i + 1 atoms encounter each 
other on nearest neighboring sites. 

For an island with radius R and infinite step edge bar- 
rier {a = 0), and in total n single atoms on top of it, 
let us approximate the encounter dynamics by a Poisson 
process, where w„(i?) denotes the encounter rate of ex- 
actly i+1 atoms. Within the Poisson approximation this 
rate can be precisely defined as the inverse average time 
for i+l atoms to encounter each other for the first time, 
when initially n atoms are randomly distributed on top 
of the island. A simple scaling argument yields 



w„(i?) 



Y[in-k) 



.k=0 



D 



-kR^ 



(24) 



where constant. The term (a} I'kR}')^^ is propor- 

tional to the probability to find \ atoms on nearest 
neighbor sites, and the factor (txR? ja}) takes into ac- 
count that the encounter can occur everywhere on the 
island. The combinatorial factor 111=0 ("'~^) slightly 
more subtle. At first sight, one may think that one should 
include the number (^j^) of possibilities to choose any 
atoms out of the n atoms but this is not correct, since 
the accumulation of atoms does not happen "in par- 
allel" at a certain instant of time but in order: First a 
dimer forms out of n single atoms (combinatorial factor 
ri(n— 1)/2) and then some of the remaining n — k atoms 
(fc = 2, 3, . . . , i) have to attach one after another to an 
intermediate cluster of size fc before this cluster disso- 
ciates (the intermediate cluster is assumed to be much 
less mobile than single adatoms). The sequential attach- 
ment process yields an additional combinatorial factor 
111=2 Clearly, the scaling argument gives only a 



I 
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FIG. 6. Scaled rate uji.a(R)(p^ I D) j X{\^^^(n-k) of the en- 
counter of i+\ atoms out of n atoms on an island with radius 
R and infinite step edge barrier. The upper curve with slope 
(-2) refers to i = 1 and n = 2 (□), 3 (o), and 4 (A), and the 
lower curve with slope (-4) refers to i = 2 and n = 3 (V), 4 (0), 
and 5 (x). 

rough approximation for ujn{R) and a more refined treat- 
ment justifying eq. ( p^ ) is presented in Appendix 

Determination of uJn (R) for i = 1 , 2 and various n in our 
simulations confirms the behavior predicted by eq. ([24|), 
see Fig. ^. For i = l the scaling law is only valid for large 
i?>100a, because at smaller R, two atoms typically en- 
counter each other before the delta-functions character- 
izing the initial occupancy smear out to a uniform distri- 
bution (for larger i this effect becomes less important). 
Moreover, we find Kq — 0.087 for i = 1 and = 0.53 for 
i — 1, i.e. the coefficient is constant for fixed i, but 
changes strongly with i. This dependence is expected, 
since we neglected the memory effect that, when n atoms, 
2<n<i, are already close to each other, they keep close 
together for a while so that the encounter of i -1-1 atoms 
during this intermediate time becomes more likely. This 
memor y e ffect is not included in the treatment in Ap- 
pendix where after each "dissociation" of an unstable 
cluster of size k<i a. configuration is assumed to emerge, 
where a cluster of size k—1 is left and the remaining n~k 
atoms are assumed to be randomly distributed. Accord- 
ingly, Ko should increase with increasing i as it is the 
case. 

Equation ( p^ ) has been derived for an infinite step edge 
barrier. For finite step edge barriers, we have to take into 
account that a state corresponding to an island with n 
atoms on top of it has a finite lifetime Tn{R) only. This 
lifetime is defined by the average time required for the 
first of the n atoms to escape from the island (if any 
encounter processes are neglected). To a good approx- 
imation, T„(i?) is the nth fraction of the lifetime ti{R) 
of a single atom, Tn{R) — Ti{R)/n (this approximation 
would become exact, if the escape were a simple Poisson 
process). In the limit of large a, ti{R) is proportional 
to the characteristic time R^ /D for an atom to reach the 
boundary, while for small a an atom typically returns 
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many times to the boundary before escaping from the is- 
land. Thus, in the latter limit, the characteristic escape 
rate (inverse lifetime t^^) is approximately given by the 
product of the probability 27ri?a/7ri?^ for the atom to be 
at the boundary and the rate aD/6a^ to overcome the 
step edge barrier. Combining these results gives 



1 f a 

Tn{R) = -— [Hl—^ 

n D \ aR 



K2 



(25) 



where ki and K2 are constants. Indeed, sua exact solution 
of the corresponding diffusion problemEZI allows one to 
derive Tn{R) exactly in the continuum limit, as we have 
shown in Appendix In particular, when the escape is 
approximated by a Poisson process, one finds ki = 1 and 
K2 = 1/2 after proper renormalization and taking into ac- 
count the lattice corrections (see Appendix ^). Direct 
determination of ti{R) in our simulations confirms this 
result, see Fig. 0. 

Knowing r„(i?) we can calculate the probability 
Pn{R) — PniR{t)) to find exactly n atoms on top of 
the island at time t before onset of second layer nucle- 
ation. This is achieved by considering the time evolution 
of Pn{R{t)), which is described by the master equation 



^^nFR{tf [il-Snfl)Pn-l-Pn] 



-r„+i(i?(i)) r„(i?(i))- 



(26) 



with the initial condition p,i(0) = (5„_o- Note that we have 
formally introduced p-i and that 1 /t„ cx n so that the last 
term on the right hand side of (^) does not contribute 
for n = 0. As can be expected and is explicitly shown 
in Appendix the solution of eq. ( p7|) is the Poisson 
distribution 



P„(i?) = ^^exp[-n(i?)], 



(27) 



where the mean number n{R) of atoms on top of the 




FIG. 7. Lifetime ti{R) of a single atom on an island with 
radius R in units of a^/D for various edge crossing probabili- 
ties a — exp{— AEs / ksT) . The solid lines are drawn according 
to eq. (bsl) with ki = 1 and K2 — I/2. 



island before onset of nucleation is 
2tt 



n{R) 



{l + aR)- 



dxx\l + 5LxY . (28) 



Here a=K2a/Ki and if = 2A^'^T'^/'^'^'^^/K2. An explicit 
solutioi i afte r evaluating the integral in eq. ( p8| ) is given 
in cq. ( B12| ) of Appendix |^. For fixed a and F, three 
distinct R regimes can be identified from eq. (p8|): For 
ipaR <C 1 we can use {l + ax)'^ ~ 1 in (28), while for 
aR <^ 1 but ipaR 3> 1 we can use (l-faa:)'^ exp{ipax) . 
For aR'3>l, we can set {l+aR)~aR in (p8|), and, since 
the integral over x is dominated by the upper bound, 
{l + ax)~ax also. We thus obtain 

-r-V(^2)^4^ i?/a«F-2/(^2)^-i 

n(i?) - <( F-ia-ii?3 , T-^i^'^'^a-^^R/a^a-^ (29) 
r-^R^, a-^^R/a 

The two regimes for large R correspond to a quasi- 
stationary situation (dpn/dt = in eq. (p6|)), where 
Pn{R) from eq. (27) equals the stationary distribution 
for R = R{t) with n{R) — TrFR^Ti{R). In these regimes 
the same result (E9) can be obtained also by integrat- 
ing pf from eq. ( fiX^ over the island area. In fact, we 
used this connection to renormalize the constants ni and 
K2 in eq. (|25|), see Appendix The small R regime 
in eq. ( p9| ) corresponds to a non-stationary situation, 
where p„ in general depends on the function R{t') at 
all times < t' < t and not only on its value R{t) at 
time t' — t. This fact, however, which also concerns the 
crossover value Rx ~F~^/^'+^)a~^ to the non-stationary 
regime, is of minor importance here, since we consider the 
generic growth law (^3|) throughout the paper. We thus 
can use R and t interchangeably. Note that the crossover 
from the non-stationary to the quasi-stationary situation 
occurs when Ti(i?x) ~ Ai(_Rx), that means in the non- 
stationary small R regime the changes in the radius occur 
on a faster scale than the escape of an atom from the is- 
land, At{R) <C Ti(i?), while in the two quasi-stationary 
large R regimes At{R) :$> ti{R). 

Let us now return to the different scenarios discussed in 
the introductory part of this Section. When n{R) > i-l-1, 
nucleation of a stable cluster can take place at any instant 
of time. The number of nucleations in At{R) that result 
from states with exactly n atoms on top of the island is 
proportional to uJn{R)At{R). The total number n„^^{R) 
is the weighted sum of cj„(i?)At(i?) over n, i.e. we find 
"nuc(^) =J2^=i+iPn{R)^n{R)'^t{li) (we are allowed to 
extend the sum up to infinity due to the sharp decrease 
of the Poisson distribution for n:^n(R)). With eq. (122" 
we thus obtain for the mean-field nucleation rate 



E 

n— frfl 

D 



Pn{R)uJn{R) 



n{R) 



t:R' 



(30) 
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Equation ( pCj ) can be interpreted as resulting from a lo- 
cal nucleation rate oc Dp^^ = D[n{R) / it B?\^^ integrated 
over the island area (factor ttR^). Compared to the 
TDT approach the radial variation of the diffusion profile 
Pi ~ Piif) is neglected in the stochastic description, aa 
that Sl(i?) from eq. ( [ll| ) may be preferred ov er cq. (^).C3 
However, as will be discussed further in Sec. IV D| below, 
for large a one should use the non-stationary solution of 



111) corresponding 



eqs. for calculating ri(i?) from 

to the small R regime of n{R) in eq. (2£ 

More important, eq. ( |30| ) (or dill)) can be used only 
if n{Rc) > i + 1 in the relevant time interval At{Rc) 
at the onset of second layer nucleation. The stochas- 
tic description allows us to treat also the fluctuation 
dominated case, where n(-Rc) <^ i + In this situa- 
tion i+1 adatoms have to be deposited and to encounter 
each other on the island. We can restrict our consid- 
eration to the deposition of exactly i + 1 atoms, since 
for n{Rc) ^ i + 1, fluctuations corresponding to more 
than i + 1 atoms on the island occur with a probability 
Y.n=mPr^{R) < exp(l)p,+i(i?)n(i?)/(z-f 2) « p.+i(i?). 
If an atom is deposited on the island already containing 
i atoms, we view this as the start of a nucleation trial. 
The number n-t^{R) of nucleation trials in time At{R) is 
ntr(-R) =TrFR^At{R)pi{R). For a trial to be successful, 
the i + 1 atoms on the island right after its start have to 
encounter each other before any of the atoms escapes by 
passing the step-edge barrier. The probability p„„^ (i?) for 
this to happen is 



p_(i?) = l-exp[-w,+i(i?)T^+i(i?)] 



(31) 



Accordingly, the total number n„„^(i?) of nucleation 
events in time At{R) is now rt„uc(i?) = nt,. (i?)penc(^), and 
using eq. ( p3|) we obtain for the fluctuation-dominated 
nucleation rate 



0,,(i?) = 7rFi?2p^(i?)p_(i?) 
= ttFR 



(32) 



■'^e-«(«)(l~exp[-c.^i(i?)r^i(i?)]) 



We note that in both formulae ( p0| , p2[) the only parame- 
ter not known a priori is the coefficient Kc, which has to 
be taken from simple simulations of the encounter pro- 
cess (see Fig. ||and the discussion above). Hence they do 
not require more input parameters than the expression 
( prj ) resulting from the TDT approach. 

It remains to clarify, when the mean-field or the fluc- 
tuation dominated situation occurs, i.e. when Q^{{R) 
or fln{R) has to be used as second layer nucleation 
rate. The answer to this question can be found by self- 
consistency requirements: Suppose first that the fluc- 
tuation dominated case takes place. Then, using (32), 
one can calculate the critical radius Rc and check if the 
condition n{Rc) <C i + 1 is fulfilled. In addition the 
condition LJi^i{Rc)At{Rc) > 1 should be fulfilled too, 
since the encounter of i -I- 1 atoms in the characteristic 
time Wi+i(i?c)~^ should happen before Rc changes. If 



these necessary conditions for the fluctuation-dominated 
case are obeyed, then the mean-field situation is ruled 
out. This conclusion can be drawn, since n{R) is 
monotonously increasing with R, which implies that Rc 
following from flfi(R) is always smaller than Rc result- 
ing from ri„f(i?). Hence, when the fiuctuations are likely 
enough to initiate second layer nucleation, they lead to 
the formation of stable clusters at an earlier time tc than 
that expected from the mean-field approach. 

We will now show that the fluctuation-dominated case 
occurs for j = 1, 2. The detailed analysis is a bit technical 
and the reader, who is interested in the main findings 
only, may skip the discussion of the various regimes I-IV 
in the following subsection and proceed with the sum- 
mary of the results given right after this discussion. 



B. Small critical nuclei {i — 1,2) 

Using fln{R) from eq. (^) we can determine the criti- 
cal radius Rc (or, more precisely, R'c) by calculating fo(t) 
as in the TDT approach (see eq. (p^)). However, for dis- 
cussing the scaling of Rc with F and a, it is easier to 
obtain Rc from the condition 



nn{Rc)At{Rc) ~ 1 



(33) 



which expresses the fact that the probability of sec- 
ond layer nucleation in At{Rc) becomes of the order of 
one. Since we consider the fluctuation-dominated case for 
small critical nuclei here (z = 1, 2), we assume n{Rc) ^ 1 
and thus set exp[— n(J?c)] — 1, when inserting f2fi(i?c) 
from eq. (^) into eq. (^). 

Four different regimes are then predicted by eq. (|33|): 

Regime I: In the limit a— +0 we have n(i?)~r^*/(*+^^i?^ 
and Ti+i (X). Hence we obtain from eqs. ( p2| , ^2p^ ) 

Fi?2r-''V(^+2)i?4^j^-lp-i/(^2)^2 _ ^4(^l)j._,(,+i)/(,^2) _ 

const., i.e. 



Rc 



p/[4(i+2)] 



(34) 



From ( |34| ) follows n{Rc) ~ const., which means that the 
assumption of a fluctuation-dominated situation is not 
necessarily justifled. In fact, eq. ( ^4| ) appears here as the 
result of a rather lengthy calculation, but in the limit 
a — ^ 0, the same scaling behavior (|3^ ) can be obtained 
very simply by calculating the average time needed for 
the deposition of z -1-1 atoms (see ref. |l6|). Hence, de- 
spite n{Rc)~i + l, eq. (34) gives the correct scaling be- 
havior. However, eq. (34) predicts uJi+i{Rc)At{Rc) ~ 

i?i?-2^j^-lr_V(^2)^2 _ p-(i^^)/[2(^2)] gi^^g p ^ 

ZJ/Fa"* 3> 1 ,Ea the inequality ujiJ^i{R)At{Rc)'>l becomes 
violated for i > 3. For i > 3 therefore, the condition 
u}i+i{Rc)At{Rc) ^ 1 should be used for calculating Rc, 
and because this yields n{Rc) > i + 1, one may alterna- 
tively us e ^l„,i( Rc)At(Rc) ~ 1 as the determining relation 
(see Sec.|IVD|). 
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Regime II: With increasing a, for i < 2, either 
the non-stationarity condition ti{Rc) ^ At{Rc) 
(n(i?c) r-*/(^2)^4 ^jjg condition 

Wi+i(i?c)Ti+i(i?c)^l (Penc — 1 in eq. (^) breaks down 
first. Taking Rc from eq. (|3^), the first condition im- 
pHes a ^ r^(*^^/[^(*+^^l, while the second imphes a <C 
p-j(2j-i)/[4(j+2)] ^ Since the first condition is more restric- 
tive for i<2, regime I ceases to be valid when a becomes 
larger than r^(*+*^/['*('+^)l and the quasi-stationarity sit- 
uation is reached. In eq. (|32) we now have to take 
n{Rc) = 7rFi?Vi(i?) - T-^oT^R^ (see eq. (Mj) and it 
follows nn{Rc)At{Rc) - a-T-*(^3)/(*f2)^34T^ ^ const., 
i.e. 



R ^ Q,V(3i44)pi(t+3)/[(*f2)(3i+4)] 



(35) 



Since n{Rc) - (r-(^)/4(^2)^-i)4/(3.+4) ^ ^ 
dition for a fluctuation-dominated situation is ful- 
flUed, and since At{Rc) ^ ti{Rc) ~ Ti+i(i?c) and 
uJi+i{Rc)Ti+i{Rc)^l the condition LJi+i{Rc)At{Rc)>l 
is obeyed too. 

Regime III: By further increasing a we obtain 
cc;,+i(i?c)r,+i(i?,)«l for « » r"'(»+-3)(2'-i)/[2(^+2)(»'-l^2)]^ 
Hence we now have to use p^^e — Wi+i(i?c)Ti-(-i(i?c) when 
inserting eq. (p2f) into eq. (p3|) and find 



(36) 



The condition n{Rf,) requires a* ^ <C 

p-(i=^-K-5)/(^2) ig fulfilled for i = 1. For i = 2, it is 
valid for a <^ r^^^"* a/Rc. The second requirement 
a;i+i(i?,)At(i?e)>l gives q,'-i < r-(*'+2''^*-5)/[('+i)(^2)] 

and again is obeyed for i — 1 and valid for i = 2 as long 
as a<r-i/4_a/i?^. 

Regime IV: In this last regime a becomes larger than 
a/ Rc, that means eq. ( |3^ ) predicts the regime to occur 
for a>r-*/[2(^2)]^ Taking n{Rc)^T-^Ri from eq. (|2|) 
and Wi+i(i?c)ri+i(i?c) i^fT^^'"^'' from cqs. (||,|5|), we 
find 



R ^ pV[2(i+2)] 



(37) 



We used Wi+i(i?c)r,+i(i?c) < 1 {p.^.{Rc) < 1) to de- 
rive (^7|), which for i = 2 is valid and for i = 1 is 
obeyed when taking into account the prefactors (for i~l 
uji+i{Rc)Ti+i{Rc) — K0K2). Moreover, eq. ( |37| ) gives 
n{Rc) ^ r'-^-^) / i^"^) ^ which is much smaller than one for 
i — For i = 2, a decision on whether the fiuctuation- 
dominated or the mean-field situation occurs would re- 
quire a closer inspection of the prefactors. However, since 
for a^a/Rc one finds the sa me scaling (^^ in the mean- 
field situation (see Sec. iVD), eq. ( p7| ) is valid in any case. 
The second condition uJi+i(Rc)AtXRc) ^ 1 is fulfilled for 
« = 1, and for i — 2 the situation again depends on the 
prefactors. 

In summary we have found that the second layer nucle- 
ation for i = 1, 2 occurs due to various mechanisms in four 



distinct regimes I- IV: In regime I (a < r~(*+*)/[^(^^)l), 
the nucleation takes place once i + 1 have been de- 
posited on the island, in regime II (r^(*+*)/['*(*+^2)l 
P_,(i4^)(2i-i)/[2(»f2)(»24^2)]) iQgg j^^Qjjjg becomes im- 
portant and the nucleation takes place once the probabil- 
ity for finding i + 1 atoms on the island at some time in- 
stant in At{R) becomes of the order of one, in regime HI 
(^p-i(i+3)(2i-i)/[2(i+2)(i=^+i+2)] <^Q,^r-'/P(^2)]^ ^Yie proba- 
bility for the encounter of i+1 atoms during a nucle- 
ation trial has to be taken into account in addition to the 
probability for the occurrence of i+1 atoms, and in regime 
IV (a 3> r^'/[^(*+^)l) both the occurrence and encounter 
probability matter but these probabilities no longer de- 
pend on the step edge barrier. For convenient reference, 
we provide the exponents 7 and n defined in eq. (|l5| ) and 
their corresponding ranges of validity in Table |. When 
comparing the scaling in the fluctuation-dominated situ- 
ation with that predicted by eqs. (|l^,^^ of the TDT ap- 
proach it is remarkable that the same behavior is found 
in regime HI and IV. We believe this to be caused by 
the fortunate circumstance that local nucleation rates of 
form cx Dp^^ might be effectively applicable even if the 
requirements for the mean-field situation are not fulfilled 
(see however ^) . 

Moreover, we have to note that for i — 1 the lower and 
upper crossovers a2 and specifying regime HI (see ta- 
ble |) both scale as T"^/^. This is due to the fact that 
when uJi^i{Rc)Ti^i{Rc) becomes less than one, we al- 
ready obtain a^a/Rc, which is the condition for regime 
IV. Nevertheless, due to the pronounced small R correc- 
tions to eq. ( [24I ) for z = 1 (see Fig. ^ both conditions 
cJi-|_i(i?c)Ti+i(i?c)^l and a-^a/Rc can be fulfilled in a 
small transient regime HI. However, for i — 1 this is no 
longer a true scaling regime, where a simple power law 
dependence of Rc on F and a can be identified. In this 
respect, the small a regime of the TDT approach does 
not occur for i — 1, also not at larger a. 



C. Comparison with simulations for i=l,2 

Taking Vta from eq. ( ^ ) we can calculate fa{t) accord- 
ing to eq. (p^). Representative results for i = 1 are shown 
in Fig. ^ (solid lines), and the comparison with the Monte 
Carlo data yields a very good agreement. The Rc values 

TABLE I. The exponents 7 and ^ characterizing the scal- 
ing Rc^V^a'^ in the various regimes I-IV for i = l,2 



Regime 


Range^ 


7 




I 


0<a<ai 


j/4(i + 2) 





II 


Ql <Q<Q2 


i(i-f3)/(i + 2)(3i-f4) 


i/iii+l) 


III 




j(j-|-3)/(i + 2)(i + 5) 


(z + l)/(j + 5) 


IV 


as <Ca < 1 


i/2{i + 2) 






^The crossover values scale as ~ F where 5\ = ij+S) /4(j+ 
2), 52 = i{i + ?.)(2i-l) /2{i + 2){f +1 + 2), and &i = i/2{i + 2). 
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derived from fo{t) are plotted as a function of a for T val- 
ues in the range 10^ — 10^ in Fig. |[ Note that, compared 
to the results of the full island model shown in Fig. |[ the 
data cover the full a range from zero to one, since the 
restrictions imposed by island coalescence in the multi- 
island model are not present i n th e single-island model 
(see also the discussion in Sec. Ill above). Moreover the 
simpler single-island model allows one to explore the be- 
havior for larger F values in the range F = 10^—10^^ also. 
It is possible to fit the Rc cu rves over the entire range of 
a and F values (see Sec. |VB| ) but we focus on the scaling 
behavior of Rc in the following in order to demonstrate 
the various scaling regimes associated with the different 
physical mechanisms of second layer nucleation. 

Indeed, the simulated data in Fig. ^ confirm the the- 
oretical predictions. For small a ^ ai (regime I), Rc 
is independent of a, while for a 3> ai ~ F"'^/'* (regime 

II) we find Rc - F^/^ia^/^. Since Rc{ai) ~ a^^^^ at the 
crossover, the boundary line between regimes I and II has 
slope (-1/9). The correctness of the scaling of Rc with F 
in regimes I and II can be deduced from the offset of the 
curves for various F in Fig. |^ (alternatively, one can col- 
lapse the data onto a common master curve by a proper 
rescaling as it was shown in ref . [l6|) . Regime II is followed 
by the transient regime III, and the dashed border line 
separating regime III from regime II was determined nu- 
merically from the condition lo2{Rc)t2{Rc) — 1 by using 
the results for lj2{Rc) and T2{Rc) displayed in Figs. pp. 
For a > r-^/[2(^2)] (regime IV), we find i?^ ~ F^/I^C^^Tl 
independent of a, and the boundary line between regimes 
III and IV has slope (-1). 
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FIG. 8. Critical island size Rc as a function of a obtained 
from simulation of the single-island model for 8 different F val- 
ues starting from 10^ (□) and ending at lO'^^ (x). Between 
these values, P is increased by a factor of 10. The various 
regimes I-IV are indicated together with the border line with 
slope (-1/9) between regimes I and II and the border line with 
slope (-1) between regimes III and IV. The dashed border line 
separating the transient regime III from regime II was calcu- 
lated numerically from the condition uj2{Rc)t2{Rc) — 1 (see 
text). 
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FIG. 9. The various regions characterizing the mechanism 
of second layer nucleation for i — 1 in an a — F diagram. 
The thick dashed line with slope (-6) marks the onset of 
layer-by-layer growth (the circles refer to the onset of island 
coalescence obtained in the full simulation, see Fig. P). 

Figure ^ depicts the various regions characterizing the 
mechanism of second layer nucleation for z = 1 in an 
a — T diagram. Varying F and a within one of the re- 
gions results in the corresponding behavior of Rc accord- 
ing to eqs. ( |3^ , |35| , p7| ) . The border line between regions 
I and II has slope (—4/3), between regions III and IV 
slope (—6), and the dashed line marks the border line 
between regions II and III. In addition, we have drawn 
the transition line from rough multilayer to smooth layer- 
by-layer growth into the diagram. In our simulations is- 
land coalescence occurs in regime II (see Fig. where 



1/2. 



'Fi/6 thus 



i?^^F4/2iQ;i/'^. The criterion Rc~l^px 
yields a^{T)r-^T-^'^ . 

Results for Rc obtained from simulations for a critical 
nucleus of size z = 2 are shown in Fig. Again the 
results confirm the predictions of the theory. In particu- 
lar, for large F, the exponents = 1/5 in regime II and 
/i = 3/7 in regime III can be clearly identified. In con- 
trast to the behavior for i = l shown in Fig. ^, regime III 
develops into a full scaling regime. 

D. Large critical nuclei (i>3) 

Analogous to the fluctuation-dominated case treated 
in the previous subsection we can obtain the scaling of 
Rc with F and a from the condition ri„,f (i?c)At(i?c) — 
1 with n^i{R) and A<(i?c) from eqs. (|0|j2|), re- 
spectively. For critical island radii belonging to 
the two quasi-stationary large i?-regimes in eq. (|29| ) 
this gives the same behavior ([l^jl^) as predicted 
by the TDT approach. However, for large step 
edge barriers corresponding to the 
small a-regime in eq. ( p9| ) we find 

Z?(F-»/('^2)^4y.H^-2«^-lp-V(«+2)^2 

const., i.e. 



non-stationary 
{Rc)At{Rc) ^ 



^2(i4-3)p 



(38) 
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10 " 10 " 10" 
a=exp[-AE5/kBT] 

FIG. 10. Critical island size Rc for i = 2 obtained from the 
single-island model for 5 different F values between 10^ (□) 
and 10^ (o). Between these values, F is increased by a fac- 
tor of 10 as in Fig. 5. The scaling regimes I-IV are indicated 
together with the respective border lines with slope (-1/5) be- 
tween regimes I and II, slope (-1/3) between regime II and III 
and slope (-1) between regimes III and IV. 

With increasing a this scaling breaks down when Rc en- 
ters the quasi-stationary regime in eq. (^9|) that means 

for a > r-2/(^2) jj-l ^ p-(i2_^ifl0)/[2(^f2)(^3)] _ 

For j > 3 we thus have in total three distinct regimes 
I-III with different mechanisms for second layer nucle- 
ation: In regime I {a < r-(''^^io)/[2(^2)(^)]) ^^le nu- 
cleation takes place once the island radius R has grown 
large enough so that the encounter of i -I- 1 atoms out of 
typically n{R) > i+1 atoms happens in a time compara- 
ble to At{R), in regime II n{R) becomes dependent on a, 
while in regime III, for large a^a/Rc, ti{R) no longer 
depends on the step edge barrier and fi[R) becomes inde- 
pendent of a again. The overall behavior characterized 
by the scaling exponents 7 and /i is summarized in Ta- 
ble m Computer simulations for i = 3 are in accordance 
with these theoretical predictions, see Fig. |l^. The pre- 
dicted scaling Rc ~ a^/^ in regime II is not yet fully 
developed for the F-values in the range 10^ — 10^ but it 
can be expected to become more clearly visible for larger 
F. However, we could not obtain reliable simulation re- 
sults for larger F values, since the amount of CPU time 
for determining the onset of second layer nucleation be- 
comes tremendous due to the increasing number of atoms 
contributing to the nucleation event. 



E. Influence of metastable clusters 

To demonstrate how the presence of metastable nuclei 
may be in cluded into the general procedure presented in 
Sec. IV A , we consider, as in ref. |l^, the simplest case of 
second layer nucleation of a trimer (i = 2), when a dimer 
is metastable with characteristic dissociation time t^^^. 
For 1 = 2 we have to deal with the fluctuation-dominated 
situation. We note in passing that this can be true even 
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FIG. 11. Critical island size Rc for i = 3 obtained from the 
single-island model for 4 different F values between 10* (□) 
and 10* (V). The scaling regimes I-III are indicated together 
with the respective border lines with slope (-5/17) between 
regimes I and II and slope (-1) between regimes II and III. 

for larger i when metastable clusters can form, since their 
presence tends to drive second layer nucleation into the 
fluctuation-dominated situation. 

In contrast to the discussion leading to eq. ( ^2[ ) for 
the non-interacting particle model, the formation of the 
stable trimer is not necessarily the rate limiting process. 
It is possible that the dissociation time t^-,^ becomes so 
large that the nucleation happens effectively instanta- 
neously once the dimer has formed. To decide whether 
the formation of the stable trimer or metastable dimer 
is rate limiting, we have to compare pi(i?)p<^'^(i?) with 
P2{R)pI^^{R), where p'£,XR) deno tes the encounter prob- 
ability of j atoms (in Sec. IV A no superscript (j) was 
introduced, since only j = i + 1 had to be considered). 
Hence we write 

r nFR^p2 (i?)pf„i (R) , pipfi > P2Pf„l 
nn{R)^< (39) 
[ ttFR^pi (i?)pf2 (R), PiPill < P2Pi'l 

To calculate the occupation probabilities Pn{R), n < 
3, we first need to know the modified lifetimes t^{R) 
of states with exactly n atoms on top of the island. 
Clearly, t[{R) = ti{R) with ti{R) from eq. (p5|), since 
the metastable dimer has no infiuence on the lifetime of 
a single atom. The characteristic time T2{R), however, 
will be enlarged in comparison to T2{R) from eq. (|25| ) 
and can be estimated as follows (we disregard any prefac- 
tors) : As in ref. |l^ we consider the first deposited atom as 

TABLE II. The exponents 7 and /i characterizing the scal- 
ing _Rc~r^Q^ in the various regimes I-III for j>3 



Regime 


Range^ 


7 




I 


0<a<ai 


(j-l)/2(i-h3) 





II 




i{i+3)/{i+2){i+5) 


(i+l)/(i+5) 


III 


a2<Ca< 1 


i/2{i + 2) 






^The crossover values scale as a*, ~ F 
10)/2(j-f2)(i+3), and S2=i/2{i + 2). 



where 5i = (i -|-5i-|- 
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immobile and the second deposited atom as diffusing. 
Once the second atom has been deposited it needs a 
time of order {R^ / D + Tdj^) to reach the step edge, since 
one encounter with the first deposited atom typicaUy 
takes place during one traversal of the island within time 
Tt, = E? /Dx]\ At the boundary the second atom is "re- 
flected" a typical number M~a~^ of times before leav- 
ing the island. Between all reflections, the overall elapsed 
time is of order {MRa/D + mTdiJ, where Ra/D is the 
typical time for a single atom to return to the edge and 
m~ {MRa/D)/Tt^^a/aR is the typical nuniber of times 
the second atom encounters the first atom.tZI Summing 
up all time contributions we obtain (neglecting the pref- 
actors belonging to the four individual terms) 



aR 



(40) 



Note that for t^,, t^, = R"^ / D , rf^iR) reduces to T2{R) 
from eq. ( psf ) (without prefactors). To estimate t^{R) 
we note that if the dimer state is the prevalent one, 
T^{R) ^ Ti{R), whereas, if all three atoms are likely to 
be separated, t^{R) ~ T3{R). Since r^lR) ^ ti{R), we 
find T^{R) ^ Ti{R) in either case. In the strong bar- 
rier limit a <C a/R, in particular, the first two terms 
on the right hand side of eq. (|4^) can be neglected, 
and, since ti Ra/Da ^ r^^a/aR, we can simply write 

~ Ti(l+rdia/Ttf) which agrees with the result derived in 
ref. This finding for strong step edge barriers implies 
that for Tdig/rtj 3> 1 the two atoms are effectively always 
in the dimer state and ~ T2Tdia/Ttj, while for T^i^/r^-^ <C 1 
they are effectively always separated and T2~r2. 

When inserting the modified lifetimes into eq. (^6|) and 
neglecting states with n>i-|-l = 3 (p„~0 for n>3 until 
onset of nucleation in the fluctuation-dominated situa- 
tion), we can calculate the occupation probabilities Pn{t). 
In the quasi-stationary limit {dpn/dt — but R — R{t)), 
in particular, we obtain (for 0<ji<i + l = 3) 



Pn 



Ejto 111=1 Ik 



(41) 



To calculate the encounter probabilities — 1 ^ 

exp[—u'„{R)T^{R)], n = 2,3, we furthermore need t o kn ow 
the modified encounter rates uj'„{R). From eq. (A6) in 
Appendix ^ we find a>2 — ''^i {^1- (A6) for i — 1) and 
tjg ~ W1W2JV0 (eq. ( A6) for i = 2), where wi ~ W2 ^ t^^"^ 
from eq. ( |AlD and v'2 = t~} (modification of eq. A2), i.e. 
Lo'2^1/Tt^ and ~ Tdi3 /Tj^^ . 

To discuss eq. ( |3^ ) we may now distinguish various 
cases depending on whether we have to consider (i) 
the non-stationary or quasi-stationary situation, (ii) the 
strong (a <C a/R) or weak barrier (a ^ a/R) hmit, 
(in) the formation of the metastable dimer or stable 
trimer as rate limiting, (iv) the encounter processes to 
be faster or slower than the escape process (w^r^ ^ 1 
or not for n = 2, 3), and (v) to be dominated by the 
metastable dimer state (t2 ~ t^^^ in the weak barrier limit 



and T'2 TiTdig/rt, in the strong barrier limit) or to be 
dominated by the state of separated atoms (t2 ~ T2). 
Rather than treating all these possible cases (and ana- 
lyzingtheir possible occurrence for the generic growth 
law (O) by employing self-consistency requirements) we 
only remark that the results obtained by Krug et alSJl are 
entailed in our description. In this work, certain regimes 
corresponding to the quasi-stationary case in the strong 
barrier limit are considered for both 91 ^ 93 <C 1 and 
9192 ^ 1, where we obtain pi ~ qi and P2 — 9192 from 



eq. (41). Since u!'it[ — Ti/t^, ~ a/aR ^ 1 in the strong 
barrier limit, we can always set p*J^'^~l in eq. (39). The 



following regimes are then discussed in ref. ^ with in- 
creasing Tdia: 

(i) For Tdi, ^ TtjTi we have p'-J^l ~ W3T3 ~ t^^^ti/t'^^ and 
T2 ~ Ti, i.e. p2 — Qi and pi ^ P2Piti- Accordingly, 
we obtain f2fi ~ {TTFR'^Ti)^TaijT^^^ corresponding to 
eq. (14) (regime I) in ref. ^ 

(ii) For t^Jti <C r^i, <C n,, we find pfj^ ~ 1 and ~ n, 
P2 — qf and pi ^ P2Pi^J^ as in (i), and hence 
~ {irFR^^Tf corresponding to eq. (15) (regime 
II) in ref |l^. In the following cases, where r^i. 



i.e. 



(Hi) 



comes even larger (and ri, Tj,. do not change), 
still have p'fj^ ~ 1 . 

For Tt,. < Tdi, < t.JttFR'^ti, ^ tit^^Jt,,, 



be- 
we 



i.e. 

P2 ~ 9192 ~ {■kFR^)'^tIt,,,Jt^,. The condition t^,, < 
t^JitFR'^tx is equivalent to pi^p2P*^^i — P2, and we 



thus find VLa ~ nFR^p2 ~ {-kFR^) 




corre- 



(iv) 



sponding to eq. (16) (regime III) in ref. 111 

For t^JitFR^ti <^ t^,^ finally, pi 3> P2 and the for- 
mation of the stable trimer is no longer the rate 
limiting process. From eq. ( ^9| ) we then obtain 
r2fi~7rFi?^pi ~ (7rFi?^)^ri corresponding to eq. (17) 
(regime IV) in ref. |l^. As expected, the scaling be- 
havior of rjf, in this limit reduces to the case i = 1 
(see also the discussion in Sec. |V C] ). 

It is clear that the above analysis is difficult to ex- 
tend to even more complicated situations. Moreover, due 
to the growing number of characteristic time scales, we 
found it increasingly difficult to discern pronounced scal- 
ing regimes in practice, see e.g. Fig. ^ We therefore 
prefer to treat the problem of second layer nucleation in 
the presence of metastable clusters within the more gen- 
eral framework outlined in the following Section. 



V. SECOND LAYER NUCLEATION IN 
GENERAL SITUATIONS 



In a more general approach to the problem of second 
layer nucleation we distinguish between different individ- 
ual states of the island during its growth with respect to 
the number of atoms that are on top of the island and the 
way a given number of atoms is decomposed into clusters 
of various sizes. Employing a Poisson approximation, the 
transition processes between the states exhibit no (in- 
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trinsic) memory and can be characterized by elementary 
rates. For the non-interacting particle model these ele- 
mentary transition rates are the deposition rate ttFR^, 
the rate for the attachments of a single atom to an inter- 
mediate cluster of size k, and the loss rate of adatoms. 
The latter is given by the inverse lifetime t^^ of a single 
atom (see eq. (p5|)). Dissociation rates enter the problem 
as additional parameters, when the lifetimes of interme- 
diate metastable clusters can not be neglected. The con- 
sequences of such dissociation rates will be discussed in 
Sec. VC. Firs t, how ever, we will present the general pro- 
cedure in part V A and show in part V B how the results 
of the simplified stochastic description in Sec. IV can be 



recovered. In addition to these previously derived results, 
it is also discussed how the general treatment allows one 
to gain detailed insight into the dominant microscopic 
pathways that are followed to form a stable nucleus on 
top of the island. 



A. General Procedure 

Let us introduce a common notation for the elemen- 
tary transition rates: Wp for the deposition rate, Wi for 

(n) 

the loss rate, W!^j for the attachment rate for a single 
atom to an intermediate cluster of size j if in t otal n 
atoms are present on top of the island (see eq. Al in Ap- 
pendix we formally include the case j = 1), and Wd.j 
for the dissociation rate of an unstable cluster composed 
of J < i atoms (again we do not distinguish between dif- 
ferent cluster configurations for the same cluster size, see 
also the remark in ref. |2^). According to the results de- 
rived in Sec. IV and Appendix ^ the transition rates are 



(42a) 
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n(n-l) 2D 
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(n-j) 



D 



(42c) 



7ri?2 ' 



2<j<n 



D 

Kd,\ exp 



2<j<i, (42d) 



where AE^'" = 



Eq is the dissociation energy of a 
single atom from an unstable cluster of size j < i. The 
prefactors k[^I and j contain the effective sizes of clus- 
ter perimeters on one hand (see the discussion in Ap- 
pendix ^) , and various corrections involved in the overall 
approximation scheme (Poisson approximation, cutoff ri-^ 
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FIG. 12. Various states and the corresponding transition 
rates (see eq. (^)) involved in the second layer nucleation for 
a critical nucleus of size i = 2. 

introduced below, etc.); they are considered to be inde- 
pendent of D, a and R. In principle one should also take 
into account the possibility that a subcluster composed 
of more than one atom can dissociate from an unsta- 
ble cluster. In fact, it has been argued that such cluster 
dissociations are sometimes more likely to occur than the 
dissociation of single atoms, as e.g. for dimer dissociation 
from a tetcamer on a (100) surface by a kind of "shear- 
ing mode" eH For simplicity we will take into account only 
the dissociation of single atoms here, although conceptu- 
ally the inclusion of cluster dissociation processes into 
the general treatment poses no difficulty. Also, we do 
not consider the influence of cluster mobilities. If one 
would allow for a small jump rate Dj/a^ of a cluster of 
size j > 2, the relative diffusion of a j cluster and a single 
atom would be larger by a factor l + Dj/D and accord- 

(n) 

ingly we had to multiply W!^ j 
forj>2. 



in eq. ( 42c ) by this factor 
The method is best introduced by an example. To this 



(42b) end consider Fig. 



12 that illustrates the situation for a 
critical nucleus of size i — 2. Various states of the island 
are shown, which are distinguished according to the total 
number n of atoms on top of the island, and the possible 
configurations that can be assumed for a given n. Be- 
tween the states the possible transitions are marked by 
arrows that are labeled by the corresponding rates. Note 
that the loss from a state with n single atoms is n times 
larger than the loss from the state with one atom. It is 
clear that Fig. 



12 shows only a small part of the possi- 



ble states and in principle can be extended by including 
larger numbers n. However, as will be pointed out be- 
low, these states with larger n do not contribute much 
to the onset of second layer nucleation. Moreover, we 
have not included states containing stable clusters of size 
j>i + l and transitions between different states contain- 
ing a stable nucleus of size i + These are irrelevant for 
the fraction fo{t) of covered islands at time t. 

We denote by Pn,u the probability for the island to be 
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in state (n, v), where n refers to the number of atoms on 
top of the island and v to a, specific configuration for a 
given n. A complete description of the stochastic process 
amounts to specifying the set {Pn,uit)} of state probabil- 
ities at all times t. The time evolution of the {Pn,u{t)} is 
described by the master equation 



dt 



\W{n' ,1^' '^n,v)pn'y 



(43) 



where for the rates Win, v n' , v') the appropriate 
expressions from eq. have to be substituted (see 

Fig. [ij) . Note that transitions are possible only between 
a limited number of states. In the situation considered 
here, where only single atoms can leave the island, we 
have W{n, v — > n', j/') = for \n — n'\ > 2. 

To treat the problem of second layer nucleation un- 
der generic growth conditions one has to solve the set of 
eqs. (HI) for R = R{t) with R{t) from eq. (|l|) subject to 
the initial condition p„_^ = <5n,o- To this end it is conve- 
nient to solve (43) using R as the independent variable. 
The integration of the differential equations (^ ) using 
standard solvers takes very little CPU time on ordinary 
workstations, so that results for fo{t) and Rc{a) can be 
obtained almost immediately. Numerical results are dis- 
cussed in the following. 



B. Negligible Lifetimes of Unstable Clusters 



In this subsection we consider the case AE'^ 



: that 

■ I J 

was treated extensively in Sec. IV. The fraction /o(i) 
of covered islands within our more general framework is 
given by 



(44) 



n—i+l 



where !/„ is the configuration containing a stable nucleus 
for a given n (for example, 1^3 = 3 and 1/4 — A in Fig. 12). 
In practice, states corresponding to large n contribute 
a negligible amount up to times tc, so that one needs 
to consider a finite maximum number of atoms n^, only 
{ni,—A in Fig. turned out to be sufficient). 

Figure 13 a, shows fa(t) and the probabilities p„^jy(t) (la- 
beled according to Fig. 12) as a function of the coverage 
Fa^t for a = 10^'*, and r = 10^. Also shown is the mean 
total number 



n—l y^Un 



(45) 



of atoms that are not in states possessing a stable nu- 
cleus. In accordance with the predictions of the simpli- 
fied stochastic description, this number is less than one 
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FIG. 13. Time development of the fraction of covered is- 
lands fo{t), the mean number of atoms N{t) in states con- 
taining no stable nucleus, the occupation probabilities Pn,i, 
and the currents js and J4 (see eq. ^ ) for i = 2, a = 10"* 
and r = 10®. The maximum number of atoms is n* = 4 corre- 
sponding to the digram shown in Fig. |l^. The vertical dotted 
line marks the coverage Fa^tc at the critical time tc. 

up to time tc- Accordingly, the pathway followed by the 
system to form a stable nucleus is dominated by fluc- 
tuations as discussed in Sec. IV. The important role of 



the fluctuations can even more clearly be recognized by 
looking at the state probabilities Pn,!^, v ^ Vn., and the 
currents 



h it) ^ W^^ P3,2 (t) , J4 (t) ^ W^> P4,2 it) (46) 

into the states containing a stable nucleus. As can be 
seen from Fig. [T^a, only the probabihties Pn,i(i) are sig- 
nificant, while the other state probabilities Pn,2it) and 
P4,3{t) cannot be discerned on the scale used in the fig- 
ure. On the other hand, we find that the current jsit) 
from the state (rt = 3, = 2) (which has a very small prob- 
ability P3.2{t)) contributes most to the growth of fo{t), 
see Fig. p^. The fact that j4{t) gives only a subdom- 
inant contribution to second layer nucleation, indicates 
that the incorporation of states with n > 5 will not sig- 
nificantly change the behavior of fo{t). 

The results for fo{t) compare well with the data ob- 
tained from kinetic Monte Carlo simulations, the qual- 
ity of agreement between theory and simulation being as 
good as in Fig. ^ The values of the optimal prefactors 



(4), 



and Kd,2 are given in ref. To exemplify the good 
agreement between theory and simulations, we have re- 
plotted in Fig. |l^ the critical radius Rc as a function of 
a for various F for i = 1, 2 from Figs. HJl^. The solid lines 
referring to the numerical results give an excellent fit to 
the Monte Carlo data. For i = 1, only the states with 
n < 2 in Fig. |l^ had to be included to achieve this almost 
perfect agreement. 
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FIG. 14. Comparison of critical island sizes Rc obtained 
from the rate-equation approach (solid lines) with Monte 
Carlo data for critical nuclei i — 1,2 and different values of F 
(same symbols as in Figs. p|Jlo|). 

For i > 2 we expect that a very large number n^, has to 
be chosen in order to obtain a correct description of sec- 
ond layer nucleation within the rate equation approach. 
Diagrams corresponding to that shown in Fig. 12 then 



become very complicated and not easily tractable from 
the practical point of view. It is thus helpful to introduce 
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FIG. 15. States and corresponding transition rates (see 
eqs. involved in the second layer nucleation for a crit- 

ical nucleus of size i — 3 when the lifetimes of unstable clusters 
are neglected. Only states corresponding to n < 5 atoms on 
top of the island are shown. 
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FIG. 16. Time development of the fraction of covered is- 
lands fo{t), the mean number of atoms N{t) in states con- 
taining no stable nucleus, the occupation probabilities Pn.i, 
and the currents jn{t) = uJnPn,!, n>i + l for i = 3, a = 10~^ 
and F — 10*. The maximum number of atoms included in a 
diagram of the type shown in Fig. |l^ is n* = 50. The vertical 
dotted line marks the coverage Fa tc at the critical time tc. 
At tc we find N{tc) = 10. 

the "renormalized" encounter rates ujn defined in eq. ( |24|) 
and to consider simplified diagrams as shown in Fig. |15| 
for z = 3. For a given number n>i + l of atoms on top 
of the island we have only included two states v = 1,2: 
One of these refers to a state where all n atoms are sepa- 
rated (j/= 1), and the other one to a state, where exactly 
i + l atoms form a stable nucleus, while the remaining 
n~{i + l) atoms are not bound to other atoms in the 
same layer {i' — 2). 

Plots of fo{t), p„,i(t), N{t), and i„(t) = 
Pn.i(t)uJniR{t)) for i = 3 analogous to Fig. O are shown 



in Fig. hq. As expected from the discussion in Sec. IV, 
we now had to take into account states with n up to 
n^ = 50^i{i + l)/2 = 6 before reaching the limit, where 
fo{t) as calculated from eq. (Q) did not change much by 
incorporation of states with larger n. Near tc, N{t) is sig- 
nificantly larger than i+l = 4 (at tc we find N{tc) = 10), 
and the dominant currents jnit) initiating second layer 
nucleation are those for n=15 — 203>4, see Fig. ^ 

In order to see how the preferred paths for second layer 
nucleation change with the step crossing probability a, 
we define the integrated current jn(t) up to tc by 



dtjnit)^ f dtpn,i{t)uJn{R{t)). 

Jo 



n>i + l. (47) 



This quantity equals the fraction of covered islands at 
time tc for which the stable nucleus originates from a 
state possessing exactly n adatoms. Figure |l^ shows J„ 
as a function of n for fixed F = 10® and various a. We 
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FIG. 17. Integrated nucleation currents up to time tc as a 
function of the number of particles n in the state from which 
the nucleation event took place (i = 3). 

see that for all a states with n ^ 4 dominate the onset 
of the nucleation. The number of particles npeak in the 
state where J„ has a maximum strongly increases with 
increasing a. For a = 10"'^, the second layer nucleation 
is typically initiated by ripoak — 18 adatoms on top of the 
island. 
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FIG. 18. Comparison of critical island sizes Rc obtained 
from the rate-equation approach with results from kinetic 
Monte Carlo simulations of the one-island model for i = 2. 
Dimers are metastable and dissociate with a rate PD/a^ (see 
eq. ©). 



C. Influence of Metastable Clusters 



VI. SUMMARY AND DISCUSSION 



The general procedure outlined in Sec. VA allows us 
also to describe situations, where the binding energies of 
unstable clusters of size j <i are not small compared to 
kBT. To demonstrate this we again consider the case i = 2 
and the corresponding diagram in Fig. |l^. The dimer in 
the intermediate states possessing no stable nucleus is 
now considered to be metastable, and we introduce the 
parameter 



(3 = exp(-A^f /fceT) 



(48) 



as "dissociation probability" (analogous to the step edge 
crossing probability a = exp{— AEs / kBT)) . For (3 = 1 
we recover the non-interacting particle model. From the 
outset it is clear that second layer nucleation will proceed 
faster for smaller /3, since the state probabilities P3,2(i) 
and Pi,2{t) in Fig. ^ and accordingly the currents j3(i) 
and jA{t) defined in cq. 46 will become strongly enhanced. 

Figure |l| shows Rc as a function of a for fixed r= 10 
and various /3 obtained from Monte Carlo simulations 
(open symbols). As expected, the critical radius de- 
creases with decreasing (3. In fact, for (3 = 10^^ one 
can regard the dimer as effectively being stable on the 
relevant time scale so that the changes with [3 corre- 
spond to a continuous transition from i = 2 {(3 = 1) to 
1 = 1 {(3 = 10"^). The comparison with the numerical 
solution of eqs. ^ (solid lines) yields very good agree- 
ment. To achieve this agreement, we used the same set 



of prefactors Kd,2 as in Fig. |l|.ii This highlights 



the power of the rate equation approach to treat second 
layer nucleation in general situations. 



In summary, we have presented a detailed theoreti- 
cal investigation of the nucleation on top of islands in 
epitaxial growth. In the non-interacting particle model, 
where the lifetimes of unstable clusters can be neglected, 
it was possible to tackle the problem within a simpli- 
fied stochastic description based on scaling arguments. 
An important result for the non-interacting case is that 
the nucleation for critical nuclei of size i < 2 is domi- 
nated by fluctuations, while for larger critical nuclei it 
can be treated in a mean-field type manner. The sec- 
ond layer nucleation rate for both cases was derived in 
compact form (see eqs. (32[3^)). When metastable clus- 
ters can form with appreciable lifetimes, the sim plifie d 
description can in principle be extended (see Sec. IV E), 
but becomes of limited value due to the fact that many 
elementary processes get mutually coupled both sequen- 
tially and in parallel. In such situations it is better to 
employ the more general framework outlined in Sec. ^ 
that is based on our derivation for the transition rates 
of the elementary processes. Results obtained from both 
theoretical approaches were shown to agree with Monte 
Carlo data. 

Throughout the paper, we have used the generic 
growth law ( p^ for the mean island radius, but it is 
straightforward to treat other growth laws also (as e.g. an 
exponential behavior), psthich may be realized by special 
preparation techniques.Eil Neither the general expressions 
(^,^2|) for the second layer nucleation rates in simple sit- 
uations nor the master equation ( ^3| ) depend on the spe- 
cific form of the growth law (the expressions for n, p„, 
etc. in the quasi-stationary case however get modified, see 
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the discussion in Sec. IVA). Moreover, it is straightfor- 
ward to rewrite all formulae for the case of heteroepitaxy 
by replacing the jump rate D/a? of adatoms on top of 
the islands by a modified jump rate D' /a? . 

The theoretical understanding of second layer nucle- 
ation is not only of basic importance but has numerous 
applications. One of these is the determination of the 
effective step edge barrier Ai?s for systems, where the 
more direct and simpler method via the measurement 
of adatom lifetimes by field ion microscopyO cannot be 
applied. As pointed out in ref. 16, the breakdown of 
the TDT approach in the fluctuation-dominated situa- 
tion calls for a reexamination of some experimental data 
for estimating Ai?s. In fact, such reexamination has been 
carried out recently by Krug et aLEZI with notable results: 
By reanalyzing the fraction of covered islands f{t) mea- 
sured for Ag/Ag(lll)EJ they corrected the previously re- 
ported estimate AE'^ w 0.12eV to AE'^ w 0.32eV (they 
also reported another estimate yielding IS.Es ~ 0.20eV 
based on a modified data analysis, see the comment in 
ref. lU). Krug et al. moreover studied the influence of 
step decoration by CO moleculea^ on Rc (and hence 
AEs) for Pt/Pt(lll). They found a strong increase of 
AEs with CO partial pressures, when analyzing the data 
corresponding to regime II (for of the fluctuation- 

dominated situation. Hence contamination by CO is ex- 
pected to favor multilayer growth. 

On the other hand, surfactants may promote smooth 
layer-by-layer growth. For example, the presence of 
only small amounts of Sb for growth of Ag on Ag(lll) 
were shpsm-to convert rough nmltilayer to layer-by-layer 
growthO£j It was suggestedEj that Sb reduces AEs, 
but, since it was observed-ihat Sb increases the island 
density in the first layerj£jiZl it is also possible that the 
induced layer-by-layer growth results from a decrease of 
the mean island distance. Even in the absence of surfac- 
tants, a change of the effective step edge barrier may go 
along with a shape transition of the islands with vary- 
ing temperature (see refs. |5^-|60|,[7| and the comment in 
ref. psf ), and this can induce changes in the film mor- 
phology as well. With respect to the transition from the 
fluctuation-dominated to the mean-field type situation 
with varying i predicted in this work, it would also be 
interesting to conduct proper experiments for metal epi- 
taxy on (100) surfaces, where a change from i = 1 to i = 3 
is often observed with increasing temperature. 

A further application pertaining to the design of self- 
organized nanostructures is the possibility to create pyra- 
midal ipoimds on a substrate, which are calle(jL."wedding 
cakes" .EJl£j As suggested by Michely et a/., a the size 
^top of the top terrace of the pyramid should be roughly 
given by ri(Ltop) ~ E, where $7 is the second layer nucle- 
ation rate. Recently, an expression for the distribution of 
Ltop has been suggested within a self-consistentpanalysis 
of a model for the dynamics of the top terrace.O In re- 
cent developments of nanostructure formation also larger 
clusters of atoms are considered as basic building blocks 
in epitaxial growth. The underlying processes seem to 



be very similar to the case of deposition of single atoms 
or simple molecules (for a recent review, see Ref. ^), so 
that it could well be that also for cluster deposition an 
effective step edge barrier has a decisive influence on the 
fllm topography. 

In light of the basic importance and the manifold appli- 
cations, there is certainly need for further improvement 
of our understanding of second layer nucleation. Topics 
worthy of further study are in particular the influence of 
strain effects and longer-range interactions between the 
adatoms. The latter may be attributed to direct forces 
(e.g. induced dipole-dipole forces in the case of magnetic 
adsorbates), or they can be mediated by perturbations 
of the electron structure of the substrate. By extending 
the approach presented in this work, these issues may be 
tackled in the near future. 
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APPENDIX A: 

We want to calculate the characteristic rate LOn {R) for 
an encounter of z -I- 1 atoms, if initially n> atoms 
are randomly placed on top of an island with radius R 
and infinite step edge barrier (a = 0) . For this purpose 
let us consider the encounter as a sequential process as 
depicted in Fig. ^ (for i — i and n — b): First a dimer 
forms, then one of the remaining atoms attaches to the 
dimer and a trimer is created, and so on until a stable 
cluster composed of i+1 atoms has been formed. Denot- 
ing the rate for the formation of the dimer by wi, and 
the rate for the attachment of an atom to an already ex- 
isting cluster composed of k atoms ( "/c-cluster" ) by w^, 
we may write 



Wi 



n(n-l) 2D bia"^ 
2 ~^ ^ 



Wk = (n-fc) 



2<fc<! 



(Ala) 



(Alb) 



The factors bk can be viewed as the effective number of 
perimeter sites of a fc-cluster. Similarly, we may write 
for the rate of dissociation Vk of a single atom from a 
fc-cluster (in the case of negligible binding energies of un- 
stable clusters) 



Vk ^ dk^ , 



fc>2. 



(A2) 
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APPENDIX B: 
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stable 

V2 V3 
FIG. 19. Illustration of the encounter of 4 atoms out of 
n = 5 adatoms. The 4 connected atoms are supposed to form 
a stable cluster {i — 3). First a dimer forms with a rate wi 
then a trimer with a rate W2, and finally the stable quadrumer 
with a rate W3. The sequential process only leads to the for- 
mation of a stable cluster, if neither the dimer dissociates with 
the rate V2 nor the trimer dissociates with the rate U3. 



The solution of the diffusion 



problenS 



dp 
dr 



dp a 

7r + -p 

or a 



(Bla) 
(Bib) 



r=R 



with the initial condition p(r,t = 0) = l/(7ri?^) has been 
derived by Harris :CJ 



where again dk has the meaning of an effective number of 
perimeter sites. (In principle one may also take into ac- 
count the possibility that a subcluster composed of more 
than one atom can dissociate from an unstable cluster 
and other states with various intermediate unstable clus- 
ters of size 2<k<i.) 

The idea now is to renormalize the process depicted in 
Fig.|l^ by replacing it by an effective transition rate Wcs 
between the initial state composed of n isolated atoms 
and the final state containing the stable cluster. Clearly, 
such a replacement is only approximately valid. After 
the replacement, the encounter rate a;„(i?) in eq. ( |2^ ) 
can be identified with Wcs- In order to derive WcS, we 
consider a stationary situation, where the probability pi 
of the initial state is kept fixed and a constant current 
J flows between neighboring states containing a k- and 
fc+l-cluster. We thus write 



^27ri?2^ •^o(Afe) 



(B2) 



Here Ji/(.) is the Bessel function of i^th order, = 
4:{aR/a)'^/{Xl[Xl + (aR/a)^]) and Afc is the fcth root 
(Ai<A2<...) of 



(^)jo(A) = AJi(A). 



(B3) 



The solution ( B2 B2 ) describes the probability density 
for a single diffusing atom that at time i = is ran- 
domly deposited on top of a circular island with a par- 
tially reflecting boundary. The probability that the 
atom has not escaped from the island up to time r is 
^'(T) = 27r drr p{r,t), which yieldS 



J = WkPk ~ Vk+lPk+l 



l<k<i-l. 



(A3) 



where pk denotes the probability of the state containing 
a /c-cluster. The set of equations (A3) can be readily 
solved for pi yielding 



J = w,p, = wipi n ^ " -^E n 



fc=2 



Wi 



k=2 j=k 



On the other hand we have 
Eliminating J from eqs. (1A4 



A5), we obtain 



WcS 



Uk=2 Wk/Vk 
I + EL2 U]=k^j/^3 



(A4) 



(A5) 



(A6) 



For large radii R^a, it holds wj/vj ^ 1 so that we can 
neglect the s um over k in the denominator on the right 
hand side of (A6). Hence we find 



]J{n-k) 



.k=0 



D 



(A7) 



where = h [11=2 / dk 



*M ^HcfeCXpf-A^^T 



fe=i 



'R^ 



(B4) 



Note that, since ^'(0) = 1, it must hold X^fe^i — 1- 

The probability that none of n independent atoms has 
escaped from the island up to time r is ^'(t)". Accord- 
ingly, the probability (f){T)dT that the first atom leaves 
the island in the time interval [r, r+dr] is 



(r).-^.-n*(r)-^*(^) 



dr 



dr 



D 

'rT^ 



J2 cji-^ 



(B5) 



D 



X Aj'^exp ~{Xl+... + X;jj^T 



from which for the average time T„(i?) = dT(j){T)T fol- 
lows: 



T„(i?) ^ n 



R^ 
D 



^1 



ii,... j„=i 



(B6) 



It is easy to show that ji^k < Ajt <ja.k, where ^ is the 
fcth zero of Jv{.). Since jj^^^ ^ (A; + i^/2 — l/4)7r for k^v, 
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the terms in the series of ( pq ) rapidly decrease with in- 
creasing jk, k = 1, . . . ,n (note that Cj depends on Xj). 
The le adin g term can be obtained by setting Cj — Sj,i 
in eq. (B4), which amounts to a Poisson approximation 
of the escape process, 4'(t) ~ exp(— AiZ?T/i?^). Within 
this approximation we obtain 



TniR) 



(B7) 



where Ai follows from eq. (B3). In the limit of small 
aR/a <C 1 one finds X\ ~ 2aR/a, while in the limit of 
large aR/a ^ 1, \\ ~ Jq ^. Combining these two hmits 
yields the interpolation formula 



1 i?" / a 
n D \ an 



(B8) 



with Ki = l/2 and K2 = 1/jo = 0.173. 

Knowing Tn{R) we can set up the master equation for 
the probabilities Pn{t) to find exactly n atoms on top 
of the island at time t in the presence of an incoming 
flux see eq. (p6|). Introducing the generating function 
(5(z,<) = ^^pP„(t)z" we obtain from eq. ( p6| ) 



ttFWQ 



1 dQ 
Ti{R) dz 



(B9) 



where R — R{t) from eq. (p^). Transforming variables 
from t to R an d defining Q{z,R) by Q{z,R = R{t)) = 
Q{z,t), eq. (B9) gives (a = 1 here) 



dR l+aR^^ ' dz 



27r 



^2pi/(ff2) 



i?^(z-l)Q, (BIO) 



where a = K2Q;/ki and ip — 2A ^T'^/'-^^y K2 Equa- 



tion (BIO) is a semi-linear partial differential equation 
of first order that can be solved by the method of char- 
acteristics. For the initial condition Q{z, i? = 0) = 1 we 
obtain 



Q{z, R) = exp [-(1 - z)n(R)] 



(Bll) 



which for p„(i?) = p„(i?(iV) = [dl^Q{z,R)/n\]z=o yields 
the Poisson distribution (ET^) with 



n{R) 



27r 



27r 



Jo 



+ axy 
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(l + ai?)"^-! (l + aR)'^^ -1 
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p+2 

In the quasi-stationary case {dtQ 
obtains 



tp+1 



in eq. (|B9|)) one 
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normal direction 



island boundary 



FIG. 20. Sketch of the island geometry in the neighborhood 
of a boundary site I, from which two jumps in the directions 
Si and S2 lead to an escape from the island over the step edge 
barrier. 
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to which (B12) simplifies for paR'til (see the discussion 
right after eq. (H)). 

Alternatively, we can determine n^^{R) by integrating 
pf{r) over the island area, which yields 



Hence we can improve the Poisson approximation by 
renormalizing the bare coefficient K2 = l/jo 1 — 0.173 to 
K2 = l/8 (note that ki = 1/2 does not change). 

Finally, in order to obtain the constants ki and K2 in 
eq. (|2^) , one has to take into account the "lattice correc- 
tions" . Let us denote by 1 the position of a lattice site and 
by Sj the nearest-neighbor bond vectors, i.e. for a trian- 
gular lattice Sj = (cos[27rj76], sin[27rj/6])a, j = 0, . . . , 5. 
The master equation describing the diffusion of a single 
adatoms on the island reads 



dt 



5j,t)-w{l,t)] , (B15) 



where w{l,t) is the pr obab ility to find the atom at lat- 
tice site 1. Equation ( B15| ) is valid as long as 1 is not 
a boundary site. In the continuum limit we can write 

E -=0 [^(1 + ' t) - ^(1' t)] - (1/2) E -=0 i^i ■ v)'^«(i, t) + 



0{a ), which yields eq. (Bla) when D is replaced by 



-Dcont = D/A (see also the remark in ref. As a conse- 
quence, one has to substitute D by Ucont in all continuum 
equations, in particular in eq. (B14), which means that in 
eq. ( p5|) (referring to the lattice simulations) one should 
take Ki=4(l/2) = 2 and k2 = 4(1/8) = 1/2. 

The value of ki is still not correct, since we have not 
taken into account the lattice correction to the parame- 
ter a. To derive this correction we consider a lattice site 
1 at the boundary. For example, one may encounter the 
situation sketched in Fig. where Sq and Si lead to 
sites outside the island (where w{l + Sj,t) = 0) and the 
remaining nearest neighbor sites I + Sj, j = 2, . . . , 5 are 
on the island. The equation corresponding to ( p315| ) then 
reads 
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^-Si:Mi+^„t)-»(u)i 

D 



i=2 



i) + ^ E [^(1 + ^1 ' - ^(1' ^)] 



^EKl + ^.,i)-«^ai)]- (B16) 



In a di scret ization of the second boundary condition 
in eq. (Bib) on a triangular lattice one has to elimi- 
nate the outer boundary points 1 + <5o and 1 + Si vi a 
the discretized version of the "bulk equation" ( p31b| ). 
This amounts to a cancellation of the term on the left 
hand side and the second term on the right hand side 
of (B16) in the continuum limit, and the replacement 

\^adw/dr. Hence eq. (B1 6|) co rresponds to the second 
boundary condition in eq. ( Blb| ), when a is replaced by 
Qfcont = 2a/\/3 in eq. dBlbj ). 

In general, k nearest neighbor sites of a boundary site 1 
can lie outside the island (fc = 1, . . . , 4). The weights how 
often such 1 occur and the way the normal direction is ori- 
ented with respect to the nearest neighbor bond vectors 
leading to the sites outside the island depend sensitively 
on the shape of the island edge. Hence, the factor 2/V3 
is only an estimate, which gives an impression on the in- 
fluence of the lattice correction to the coefficient ki. Our 
comparison with the simulation results in Fig. |^ yields 
acont — 2a, i.e. ki = 4(1/2)^ = 1. We note that in general 
lattice corrections have always to be included in a con- 
tinuum description after the effective Schwoebel barrier 
(for the lattice) has been calculated from the microscopic 
barriers (see the comment in ref. p3| ). 
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